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NOMENCLATURE 

A 
n 

Coefficient  of  the  n—  term  of  the  summation, 
defined  with  Eq.  (B-10) 

C 

Concentration  of  tagged  material  at  z 

Co 

Concentration  of  tagged  material  in  flask  I 
at  zero  time 

Cl 

Concentration  of  tagged  material  in  flask  II 
and  capillary  tube  at  zero  time 

c» 

Concentration  of  tagged  material  in  system  at 
equilibrium 

c 

Number  of  constraints 

D 

General  diffusion  coefficient,  cm  /sec 

D' 

Approximate  diffusion  coefficient  defined  in 
Eq.  (B-7) ,  cm2/sec 

DBi 

Diffusion  coefficient  of  species  i  at  atmospheric 
pressure,  cm2/sec 

DSiJ 

Effective  diffusion  coefficient  (refers  to  diffu- 
sion in  transition  range)  of  species  i  in 
species  j ,  cm2 /sec 

DKi 

Knudsen  diffusion  coefficient  of  species  i, 
cm2/sec 

°iD 

Diffusion  coefficient  of  species  i  in  species  j 
according  to  bulk  diffusion  theory,  cm2/sec 

F 

2 

Pressure-induced  flow  term,  g-moles/cm  -sec 

fD 

First  two  terms  of  Taylor's  series  expansion 
for  the  partial  of  Re  with  respect  to  D 

fi(r,vi,t) 

Distribution  function 

2 
Mean  molecular  flow,  molecules/cm  -sec 

J 

Kn 

Knudsen  number,  l/r 

k 

Boltzmann  constant  for  a  molecule 

k. 

i 

Proportionality  constant 

vii 

<» 

L 

Effective  length  of  capillary  tube,  cm 

L' 

Half-length  of  capillary  tube,  cm 

L[f  (t)] 

Laplace  transform  of  f(t) 

L_1[f (s)] 

Inverse  Laplace  transform  of  f(s) 

I 

Mean  free  path  length,  cm 

Mi 

Molecular  weight  of  species  i,  gm/g-mole 

m. 
1 

Molecular  mass  of  species  i ,  gm 

* 

Ni 

Molecular  transport  of  species  i  in  molecules 
per  unit  time 

(N-c) 

2 
Number  of  degrees  of  freedom  associated  with  o 

ni 

Number  of  molecules  of  species  i  per  unit  volume 

P 

Pressure 

R 

Gas  constant  per  gram 

Re 

Residue 

r 

Radius  of  capillary  tube,  cm 

s 

Laplacian  variable 

T 

Absolute  temperature,  °K 

* 

Temperature  dependent  potential  energy  parameter 
for  i-j  interaction,  (e. ./k) 

Tl/2 

Half-life  of  pressure  gradient  decay,  sec 

t 

Time 

Vl 

Volume  of  tagged  reservoir,  cm 

V2 

Volume  of  untagged  reservoir,  cm 

V 

Arithmetic  mean  speed,  cm/sec 

K(S,t) 

Dimensionless  concentration  variable,  (C/CQ) 

Xci 

Calculated  value  of  X  at  time  i 

Xei 

Experimental  value  of  X  at  time  i 

viii 

X. 
la 

Approximated  value  of  X  at  time  i 

Xo, 

Equilibrium  value  of  X 

la   lb 

Mole  fractions  of  species  1  at  each  end  of  the 
diffusion  path 

z 

Distance  variable  along  tube 

a 

Dimensionless-parameter  for  volume  of  tagged 
reservoir,  (R  L/V.) 

a' 

Defined  with  Eq.  (20) 

S 

Dimensionless  parameter  for  volume  of  untagged 
reservoir,  (R2L/V2> 

iDL 

Correction  of  the  diffusion  coefficient  for 
finite  tube  length,  cm2/sec 

AXd 

Change  in  the  normalized  activity  in  flask  I 
by  diffusive  flow 

Total  change  in  the  normalized  activity  in 
flask  I 

AXT 

AXF 

Change  in  the  normalized  activity  in  flask  I 
by  pressure-induced  flow 

Y 

Euler's  constant 

Eij 

Force  constant  for  i-j  interaction 

Xn 

Eigenvalues  of  transcendental  equation 

w 

Gas  viscosity,  poise 

C 

Dimensionless  distance  variable,  (z/L) 

0 

Density 

0 

Molecular  diameter 

°D 

Standard  deviation  for  D 

.*<■) 

Variance  for  D  obtained  by  least  squares  analysis 

3D(P) 

Variance  for  D  caused  by  uncertainties  in  pressure 
measurements 

1 

Molecular  diameter  of  species  i 

1] 


Variance  of  pressure  measurements 

Standard  deviation  of  least  squares  calculation 
of  X(£,t) 

Variance  of  X(£,t)  calculation  by  least  squares 
analysis 

Effective  diameter  of  1-2  interaction 


2 

t  Dimensionless  time  variable,  (tD/L  ) 

e  ,<(>  Coefficient  variables  defined  with  Eq.  (A-22) 

2 

li  Net  flux  of  species  i  ,  gm/sec-cm 

nimi 
Q.  .  Potential  energy  parameter  for  i-j  interaction 


1.0   INTRODUCTION  AND  THE  THEORY  OF  GASEOUS  DIFFUSION 
1.1   Types  of  Gaseous  Diffusion 

A  vast  amount  of  scientific  interest  and  effort  has  been  and 
continues  to  be  focused  on  transport  phenomena — the  transport  of 
momentum,  energy,  and  mass  within  matter — for  their  common  and 
almost  universal  occurrence  necessitates  an  understanding  of  them 
to  ensure  the  success  of  myriad  undertakings  in  basic  and  applied 
science  and  engineering.  While  the  phenomena  are  quite  analagous 
in  many  respects,  their  differences  are  sufficient  to  warrant  and 
require  their  individual  study. 

The  similarities  of  the  phenomena  are  demonstrated  by  the 
conventional  definitions  of  the  transport  coefficients.   The 
coefficients  of  viscosity,  thermal  conductivity,  and  diffusion 
are  defined  as  the  negative  ratios  of  the  net  fluxes  of  momentum, 
energy,  or  mass  to  the  gradients  of  these.   The  definition  of  the 
diffusion  coefficient  in  a  macroscopic  system  in  one  dimension  in 
equation  form  is 

dn.  m. 

n.m.       dz 

2 

where  t|i     =  the  net  flux  of  species  i,  gm/sec-cm 

2 

D  =  the  diffusion  coefficient,  cm  /sec 

n.m.  =  the  concentration  of  species  i,  gm/cm 


n.  =  the  number  density  of  species  i,  molecules/cm 
m.  =  the  mass  of  a  molecule  of  species  i,  gm/molecule 
z  =  the  distance  variable  along  which  diffusion  occurs 


The  particular  topic  of  this  study  is  gaseous  diffusion  in 
the  transition  pressure  range,  one  of  the  three  broad  subdivisions 
of  gaseous  diffusion.   The  other  subdivisions  are:  (1)  Ordinary, 
or  bulk  diffusion;  and  (2)  Knudsen  diffusion.   A  brief  discussion 
of  each  subdivision  follows  with  an  outline  of  the  kinetic  theory 
of  gases  whose  concepts  are  common  to  all  three. 

1.2   Ordinary  Gaseous  Diffusion 

In  any  real  gas,  the  molecules  travel  in  all  directions  with 
a  wide  distribution  of  velocities,  and  molecular  interactions  are 
quite  complex  because  of  long  range  attractions  and  short  range 
repulsions.   However,  a  surprisingly  useful,  although  oversimpli- 
fied, model  for  the  behavior  of  a  gas  containing  n  molecules  per 
unit  volume  may  be  formulated  based  on  the  following  assumption: 

1)  The  molecules  are  rigid,  non-attracting  spheres  with 
diameter  c. 

2)  All  the  molecules  travel  with  the  same  speed;  a  reason- 
able choice  is  the  arithmetic  mean  speed;  for  a  Maxwellian 
speed  distribution,  v  =  (8kT)  '  (ttm)-  '   cm/sec. 

3)  All  the  molecules  travel  in  a  direction  parallel  to  one 
of  the  coordinate  axes. 

Because  of  the  first  assumption,  collisions  between  molecules 
are  well  defined,  and  a  quantity  known  as  the  mean  free  path,  the 
average  distance  traversed  by  a  molecule  between  collisions,  may 
be  introduced.   By  use  of  the  other  assumptions  and  the  ideal  gas 
law  (p=nkT) ,  this  quantity  may  be  shown  (36)  to  be 


II   =   kT/kjPwo2       ,       cm  (2) 

where  k,    =    (1/3)    +    (2/3)  /2 

(If  the  speed  distribution  were  assumed  Maxwellian  with  molecular 

motion  in  all  directions,  the  result  would  be  the  same  except  that 

—  8 
k.  would  be  equal  to  /2) .   For  simple  molecules  o  is  about  3x10 

cm  and  at  one  atmosphere  pressure  and  2  0  degrees  C.  the  mean  free 
path  of  a  gas  molecule  is  approximately  given  by  I  =  10   /p  cm 
where  p  is  in  atmospheres.   It  can  be  shown  (36)  that  by  the  sim- 
plified kinetic  theory  the  diffusion  coefficient  is  given  by 

Dll  =  v  H/3  =  k2(TrmkT)1/2Ua2p)"1  ,  cm2 /sec  (3) 

where  k.  =  2/(3n) 

p  =  Pm/kT,  the  density  of  the  gas. 

Although  the  mean  free  path  does  not  appear  directly  in  a  rigorous 
development,  the  solution  is  of  the  same  form.   However,  a 
rigorous  derivation  shows  that  k_  should  be  3/(8tt).   By  inserting 
this  value  and  combining  constants,  the  expression  becomes 

D1;l   =  2.628  x  10"3  T3/2(M1/2P''2)~1  ,    cm2/sec         (4) 

It  should  be  noted  that  this  expression  is  for  self-diffusion,  the 
inter-diffusion  of  particles  of  the  same  mass  and  size.   Using  the 
same  rigid  sphere  model  the  diffusion  coefficient  for  a  mixture 
of  two  species,  D12'  according  to  Hirschfelder  (36)  is 


D12  =  2.628  x  10~3  T3/2(y1+N2)1/2(2M1M2)"1/2(Po12)"1  ,  cm2/sec   (5) 


where  a-2  =  (a.+a2)/2   ,   cm 

As  would  be  expected,  although  the  rigid  sphere  model  yields 
reasonable  results,  these  results  give  only  approximate  pressure 
and  temperature  dependence  of  the  diffusion  coefficient  because 
of  the  severe  restrictions  entailed  by  the  simplicity  of  the  model. 
Rigorous  kinetic  theory  allows  prediction  of  these  dependencies 
with  considerably  greater  accuracy  by  considering  the  effects  of 
interactions  which  take  place  between  real  molecules.   For  most 
purposes  it  is  most  convenient  to  describe  these  interactions  by 
use  of  potential  energy  functions  which  describe  the  potential 
energy  of  interaction.   Two  of  these  functions  which  are  commonly 
used  (36)  are  the  Lennard- Jones  potential  for  non -polar  molecules 
and  the  Stockmayer  potential  for  polar  molecules.   While  these  two 
functions  are  idealizations  of  the  true  energy  of  interaction, 
many  of  the  properties  of  gases  (and  liquids)  have  been  calculated 
in  terms  of  them,  and  they  are  reasonably  adequate  for  a  number 
of  simple  molecules. 

The  rigorous  development  of  the  kinetic  theory  of  gases  is 
based  on  a  knowledge  of  the  distribution  function,  fi(r,v. ,t), 
which  represents  the  number  of  molecules  of  i    species  which  at 
time  t  lie  in  a  unit  volume  element  about  the  point  r  and  which 
have  velocities  within  a  unit  range  about  v. .   Under  equilibrium 
conditions,  f.(r,v. ,  t)  reduces  to  the  Maxwellian  velocity  dis- 
tribution (36).   When  the  system  is  not  at  equilibrium,  the 


distribution  functions  satisfies  the  Boltzmann  integro-dif ferential 
equation. 

For  the  usual  definition  of  the  diffusion  coefficient  to 
apply,  i.e.,  for  the  derivatives  of  the  flux  vectors  to  be  linear 
a  gaseous  system  must  be  near  equilibrium.   In  this  limit  the  dis- 
tribution function  is  nearly  Maxwellian,  and  the  Boltzmann  equation 
can  be  solved  by  perturbation  methods  of  Chapman  and  Enskog  (12) . 
Expressions  for  mass  flux  and  the  diffusion  coefficient  may  be 
obtained  from  this  solution. 

These  expressions  show  that  mass  transfer  may  be  produced 
not  only  by  a  concentration  gradient,  but  also  by  a  temperature 
gradient.   This  fact  is  not  displayed  in  the  derivation  by  the 
simple  kinetic  theory. 

Neglecting  temperature  effects,  the  first  approximation  to 
the  diffusion  coefficient  for  a  binary  mixture  by  the  Chapman- 
linskog  theory  is 

D12  =  2.628  x  10~3  T3/2(M1+M2)1/2(M1M2)~1(P  j^2fJ12T*2)_1  ,  cm2/sec 


(6) 


2 

where     Di?  =  diffusion  coefficient  in  cm  /sec, 

P  =  pressure  in  atmospheres, 
T  =  temperature  in  degrees  K, 

T12  =  kT/E12' 
M. My   =  molecular  weights  of  species  1  and  2, 

respectively,  gm 
°l?'  El?/k  =  molecular  potential  energy 

parameters  characteristic  of  1-2 


interaction  in  A  and  degrees  K, 
respectively, 
£2, 2  =  a  complicated  potential  energy  parameter 

which  indicates  the  deviation  of  any  part- 
icular molecular  model  from  the  idealized 
rigid-sphere  model. 

When  the  above  equation  is  written  for  a  single  component,  the 

coefficient  of  self-diffusion  is  given  by 

Dn  =  2.628  x  10"3  T3/2(M1/2Po^1n11T*1)"1  ,  cm2/sec  (7) 

2 

where  D..  =  coefficient  of  self -diffusion  in  cm  /sec,  and  the 

other  parameters  are  the  same  as  before. 

It  should  be  noted  that  here  the  diffusion  coefficient  is 

3/2 

proportional  to  I  '   and  inversely  proportional  to  the  pressure 

as  it  was  in  the  expression  derived  from  simple  kinetic  theory. 
However,  these  proportionalities  are  now  modified  by  the  presence 
of  El  and  T*  which  adjust  the  coefficient  for  molecular  inter- 
actions. 

The  Chapman-Enskog  theory  of  gases  has  generally  been  adopted 
and  is  considered  sufficient  for  this  study  for  use  in  the  investi- 
gation of  gaseous  systems  under  limited  conditions.   The  causes 
of  these  limitations  are  assumptions  in  this  treatment  which  re- 
strict the  applicability  of  the  results.   Since  only  binary  col- 
lisions of  molecules  are  considered,  the  theory  is  not  applicable 
at  elevated  pressures  where  the  frequency  of  collisions  of  three 
or  more  molecules  at  once  becomes  appreciable.   Because  of  the 


use  of  classical  mechanics  in  the  derivations,  the  theory  is 
inapplicable  at  low  temperatures  where  relativistic  phenomena 
become  significant  considerations.   The  Chapman-Enskog  method  for 
solution  of  the  Boltzmann  equation  provides  a  series  approximation 
to  the  distribution  function.   The  first  approximation,  which  de- 
fines the  diffusion  coefficient  conventionally  in  terms  of  the 
first  derivative  of  the  concentration  is  applicable  only  when 
this  gradient  is  small.   Higher  approximations  provide  corrections 
for  larger  gradients,  but  thus  redefine  the  diffusion  coefficient. 
Quantitatively,  the  first  approximation  is  valid  only  when  the 
relative  changes  in  the  concentrations  in  the  distance  of  a  mean 
free  path  are  small  compared  with  unity.   A  further  assumption  in 
development  of  the  theory  is  that  the  dimensions  of  the  containing 
vessel  are  large  compared  to  the  mean  free  path,  a  condition  which 
is  not  met  by  rarefied  gases  in  moderately  small  containers.   Since 
this  study  is  concerned  primarily  with  rarefied  gases,  these  two 
last  restrictions  necessitate  the  investigation  of  other  models 
for  the  treatment  of  low  pressure  phenomena.   Finally,  the  Chapman- 
Enskog  theory  of  gases  applies  only  to  monatomic  gases,  strictly 
speaking;  however  the  diffusion  process  is  not  appreciably  dif- 
ferent for  polyatomic  molecules,  provided  the  deviation  from  a 
spherical  shape  is  small. 

1.3   Knudsen  Diffusion 

As  the  pressure  of  a  gaseous  system  is  decreased,  the  fre- 
quency of  collisions  of  molecules  with  each  other  decreases  and 


the  mean  free  path  of  the  molecules,  j, ,  increases.   The  proba- 
bility of  a  molecule  reaching  a  wall  of  the  containment  vessel 
thus  increases  and  the  effects  of  intermolecular  collisions 
eventually  become  negligible  when  compared  with  those  of  molecule- 
wall  collisions.   This  situation  occurs  when  the  Knudsen  number, 
K  ,  which  equals  i/r  for  a  tube  of  radius  r,  is  greater  than  10. 
Under  this  condition,  there  is  little  mechanism  for  the  establish- 
ment of  local  equilibrium  and  free-molecule  or  Knudsen  flow  occurs. 
This  flow  is  analogous  to  a  process  of  diffusion  and  occurs  for 
each  constituent  along  the  gradient  of  its  partial  pressure  or 
its  concentration  gradient. 

If  the  pressure  of  the  system  is  indeed  low  enough  to  make 
the  effects  of  intermolecular  collisions  negligible,  the  charact- 
eristics of  this  flow  may  be  determined  by  consideration  of  molecule- 
wall  interactions  only.   Knudsen  (47)  experimentally  determined 
that  the  reflection  of  gas  molecules  from  a  glass  wall  is  almost 
100%  diffuse  and  follows  the  cosine  reflection  law.   He  assumed 
that  this  was  caused  by  irregular  roughness  of  the  surface  on  the 
atomic  scale.   However,  the  work  of  Langmuir  (48,49)  and  Volmer 
and  Estermann  (76)  demonstrated  that  gas  molecules  are  not  re- 
flected by  a  simple  impact  and  rebound  mechanism.   Rather,  gas 
molecules  are  absorbed  on  impact  and  after  a  finite  time  are 
evaporated.   Thus  every  surface,  no  matter  how  smooth,  will  in 
theory  reflect  gas  molecules  according  to  a  cosine  law. 

If  the  actual  irregular  container  wall  surface  is  replaced 
by  a  smooth  mathematical  surface  and  if  it  is  assumed  that  the 
adsorbed  film  is  in  equilibrium  with  the  gaseous  phase,  the 


number  of  molecules  desorbed  may  be  obtained  by  proceeding  as  if 
the  desorbed  molecules  originated  in  an  equilibrium  gas  layer  on 
the  other  side  of  the  surface.   The  number  of  molecules  desorbed 
from  a  unit  area  of  the  wall  per  unit  time  in  a  particular  direc- 
tion is  then  equal  to  the  number  crossing  the  unit  area  of  the  mathe- 
matical surface  in  that  direction.   This  quantity  is  obtained  by 
integrating  over  all  the  speeds  of  the  molecules  within  some 
solid  angle  about  the  direction  of  interest.   Using  this  expres- 
sion for  the  distribution  of  the  molecules  moving  in  a  particular 
direction,  the  net  transport  through  a  cross  section  of  the  con- 
tainer can  be  calculated.   For  a  long  tube,  Present  (60)  shows 
that  this  quantity  is 

N*  =  -  (2nr3  v/3)  dn/dz  (8) 

or  converting  to  mean  flow  per  unit  area  per  unit  time 

J  =  -  (2v  r/3)  dn/dz   .  (9) 

Equation  (9)  shows  that  Knudsen  flow  has  the  characteristics 
of  diffusion — it  is  proportional  to  the  negative  concentration 
gradient — and  the  diffusion  coefficient  is  seen  to  be 

D„  =  2  v  r/3  ,  cm2/sec   .  (10) 

K 

This  is  quite  similar  to  the  bulk  diffusion  coefficient,  (v  1/3), 
and  as  would  be  expected  the  mean  free  path  length ,  4 ,  has  been 
replaced  by  the  diameter  of  the  confining  tube. 
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It  is  noted  that  in  Knudsen  flow,  the  diffusion  coefficient  DR  is 
proportional  to  T  '  ,  to  the  pore  diameter,  and  is  independent  of 
the  total  pressure.  Using  an  assumed  Maxwellian  velocity  distri- 
bution the  result  becomes 

DR  =  (2r/3)  (8kT/Trm)1/2 


9700  rCT/M)1/2  cm2/sec   ,  (11) 


2 

Where  D„  =  coefficient  of  Knudsen  diffusion  in  cm  /sec, 

K 

r  =  pore  radius  in  cm, 

M  =  molecular  weight  of  the  species,  gm 

T  =  absolute  temperature  in   K. 

For  a  pore  diameter  of  1.026mm  at  20°C  and  at  a  pressure  for  which 

2  2 

K   =10,  D„  equals  1284  cm  /sec  for  carbon  dioxide  and  2127  cm  /sec 
n         K  ^ 

for  methane.   For  the  diffusion  of  C02  in  N2  at  25°C  in  a  pore  of 
radius  0.05  cm,  D„  and  D. ,  become  the  same  at  a  pressure  of  89.2 

K         11 

microns.   For  diffusion  at  one  atmosphere  pressure  the  coeffi- 

o 

cients  are  equal  for  a  pore  size  of  588  A. 

If  Knudsen  flow  of  a  binary  gas  mixture  occurs,  the  component 
gases  will  diffuse  along  the  tube  independently  and  Eq.  (8)  can 
be  applied  to  each  gas  separately  if  n  is  assumed  to  refer  to 

one  species  of  molecules  in  the  mixture.   For  steady  flow,  n  must 

dn 
be  constant  along  the  tube;  therefore,  the  density  gradient  g^-  is 

uniform  and  can  be  replaced  by  An/L  where  An  is  the  difference  in 

density  between  the  ends  and  L  is  the  tube  length.   Since  An  = 

Ap/kT  from  Eq.  (8),  the  flow  rate  in  moles  is 
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2 
F(moles)  =  (D„nr  )  AP/LRT   ,  moles/sec  (12) 

where  F  applies  to  the  component  with  partial  pressure  drop  Ap, 
and  R   is  the  gas  constant  per  gram.   It  is  noted  that  the  flux 
of  a  particular  component  is  inversely  proportional  to  its  mole- 
cular weight.   These  formulae  for  Knudsen  flow  have  a  sound  the- 
oretical basis  and  have  been  verified  by  experimental  work. 

1.4   Transition  Range  Diffusion 

As  has  been  discussed,  the  mechanism  of  Knudsen  gaseous 
diffusion  is  well  understood  and  rigorously  explained  as  is 
Poiseuille  or  bulk  diffusion.   The  important  range  of  transition 
from  Knudsen  to  bulk  diffusion  encompasses  a  thousandfold  range 
of  mean  free  paths.   Typical  chemical  process  conditions  fall  in 
this  transition  range  when  common  pelleted  catalysts  with  a  pore 
diameter  of  a  few  millimicrons  are  utilized.  Increased  interest 
in  ballistics  problems  arising  from  flight  in  the  upper  atmosphere 
has  accentuated  the  need  for  further  investigation  of  this  inter- 
mediate region.   The  small  pore  diameter  of  the  catalyst  pellets 
produces  a  large  value  for  the  Knudsen  number,  as  do  the  long  mean 
free  path  lengths  of  the  upper  atomsphere.   In  either  case,  Knudsen 
effects  and  transitional  behavior  result. 

The  net  transport  through  a  cross  section  of  a  long  tube  under 
conditions  of  Knudsen  flow  is  given  by  Eq .  (8).   Pollard  and 
Present  (59)  state  that  in  the  high  pressure  limiting  form  at 
uniform  total  pressure  the  analagous  expression  for  self  diffusion 
is 

N*  =  Ur2Du/kT)  dPj/dz.  (13) 


f/2 

(r/1)    = 

cosede 

. 

0 
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In  the   absence  of  a  total  pressure  gradient  and  the  associated 
molecular  drift  velocity,  the  diffusive  transports  which  may 
occur  are  amenable  to  kinetic  theory  treatment  so  that  it  should 
be  possible  to  discuss  theoretically  the  transition  from  Eq.  (8) 
to  Eq.  (13)  .   Pollard  and  Present  show  that  for  a  tube  of  infinite 
length  considering  both  bulk  and  Knudsen  self-diffusion 

N*  =  -Ur2vV3)  [l-31/(8r)  +  (61/  (itr)  )  Q  (r/U  Idn^/dz      (14) 

where 

rn/2 

-2rcosecsc*/Jl  _„2.  .  2  ...       ,1C, 
e  cos  ijisin  \Jidi^  .     (la) 

0 

Although  the  evaluation  of  Q(r/i)  often  requires  numerical  tech- 
niques, an  expression  which  can  be  used  for  small  values  of  r/J. 
is 

Q(r/i)    =   if/16    -   Trr/(6t)    +   Ttr2/(3fl2) 

-    (it/3)[1.2264    -    (3/4)ln(2Yr/Jl)  ]      r3/l3   +    ...  (16) 

where  y   is  Euler's  constant. 

From  Eq.  (14)  the  diffusion  coefficient  is  given  as 

D    =  (v  1/3)  [l-3*/(8r)  +  (6t/(nr))Q(r/t)]  .    (17) 
11 

Pollard  and  Present  further  state  that  for  a  finite  tube  of  half 
length  L1  the  diffusion  coefficient  must  be  reduced  from  that  for 
an  infinite  tube  by  a  factor 
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ADJ  =  (v  r/2)  [e~L'/l1  +  (L  '/l)Ei  (-L  70  ]  (r/L ' )   .         (18) 
L 

This  correction  term  is  negligible  for  mean  free  paths  a  few  times 
shorter  than  the  half  length  of  the  tube. 

Eq.  (17)  agrees  very  closely  with  an  earlier  expression  ad- 
vanced by  Bosanquet  (6) 

D_1   =  Dl1  +  D7|  (19) 

ell     Rl     U 

which  states  that  the  resistance  to  transport  is  the  sum  of  resis- 
tance caused  by  wall  collisions  and  intermolecular  collisions.   It 
is  derived  by  considering  the  diffusion  process  in  the  tube  as  a 
random  walk  process  in  which  the  successive  steps  of  the  individual 
molecules  are  terminated  either  by  collisions  with  other  molecules 
or  with  the  tube  wall.   The  mean  step  size  is  then  related  to  the 
mean  free  paths  for  wall  collisions  and  gas  collisions  by  taking 
the  total  collision  frequency  to  be  v/l.      The  additive  resistance 
law  results  then  from  knowledge  that  the  frequencies  of  the  two 
types  of  collisions  are  additive  and  that  D   and  D..  can  be  taken 
proportional  to  v  times  the  mean  step  size.   This  relationship  is 
probably  the  most  widely  accepted  for  the  diffusion  coefficient 
in  the  transition  pressure  range. 

Scott  and  Dullien  (6  7)  support  Bosanquet' s  expression  through 
a  derivation  based  on  the  Stefan-Maxwell  momentum  balance  method 
(60)  in  a  binary  system.   A  momentum  balance  on  the  molecules  of 
species  1  is  made  in  a  differential  volume  element  by  equating  the 
net  flux  of  momentum  carried  in  and  out  of  the  volume  element. 
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From 

this  balance  an  expression  for  the  net  mole 

cular  transport 

of  1 

is  derived 
N*  =  (P/RTL)  (D12/a')ln{(l-a,ylb  +  D^/C 

■     (1-°'yia  +  D12/DK1)_1} 

fa* 

(20) 

where  a'  =  1  +  N*/N*  , 

y,   and  y. .  refer  to  the  mole  fraction  of  1 
'la     Jlb 

at 

the  ends 

of 

the  diffusion  path. 

By  equating  N*  as  derived  from  momentum  with 

N* 

as  < 

jiven  by 

a 

usua 

1  simple  diffusion  coefficient 

Nl  "  De12P/RTL(yia-yib> 

(21) 

the 

effective  diffusion  coefficient  is  shown 

D     =  D,,/xT„ 
e,„     12'  LM 

to 

be 

(22) 

where  X.  ,  is  the  logarithmic  mean  average  of 

Xl 

and 

x«  , 

XLM  "  <X2-X1)lnX1/X2  =  a'(yla-ylb)ln{(l- 

r«',rx.+D12/Htl) 

•<1-»'yib  +  Di2%,) 

(23) 

The 

definitions  of  X,  and  X2  are 

Xl  =  X  "  -'^la  +  D12/^l 

* 

and 

X2  =  l    ~    a'ylb  +  D12/DK1 

■ 
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For  self  diffusion,  a'  is  zero,  and  the  expression  of  D    re- 

ell 

duces  to  that  of  Bosanquet. 

Rothfeld  (66)  presents  a  derivation  similar  to  that  of 
Scott  and  Dullien.   Evans,  Watson,  and  Mason,  (22)  also  support 
Bosanquet' s  relation  by  a  rather  unusual  treatment.   Theirs  is 
a  model  in  which  a  porous  solid  diffusion  medium  (such  as  a 
catalyst  pellet)  is  visualized  as  a  collection  of  spherical 
particles  which  are  very  large  in  size  as  compared  to  a  gas  mole- 
cule.  The  system  is  then  treated  as  a  "dusty"  gas  in  which  the 
dust  particles  are  held  stationary.   By  formally  treating  the 
dust  oarticles  as  giant  "molecules"  which  are  merely  one  com- 
ponent of  a  multicomponent  mixture,  the  problem  becomes  one  of 
classical  kinetic  theory  of  multicomponent  mixtures. 

Wheeler  (82)  proposes  an  exponential  function  to  define  D 

D     =  D..  (1  -  e"Dll/DKl)    .  (24) 

ell     1X 

However,  his  reasoning  is  primarily  intuitive  and  the  relation- 
ship admittedly  semi-empirical.   This  is  the  only  expression 
founa  in  the  literature  which  deviates  significantly  from 
Bosanquet's  relationship.   For  CO,  and  CH.  systems  at  20°C  this 
effective  diffusion  coefficient  is  shown  in  Fig.  1  as  a  function 
of  pressure  for  diffusion  through  a  1.026  mm  I.D.  tube.   It  is 
noted  that  these  combinations  rules  disagree  in  the  range  of  three 
to  three  hundred  microns  of  mercury.   This,  in  essence,  comprises 
the  transition  range. 
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Figure  1.   Theoretical  effective  self  diffusion  coeffi- 
cients for  CO2  and  CH4  at  20°C  with  a  0.1026  cm  I.D. 
tube  (from  Eqs.  (7,  11,  19,  and  24)). 
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2.0   EXPERIMENTAL  THEORY 


2.1   Introduction 


The  experimental  arrangements  used  to  measure  gaseous  dif- 
fusion coefficients  all  have  in  common  the  establishment  of  a 
concentration  gradient  for  the  component  of  interest.   These 
arrangements,  however,  are  quite  distinct  in  other  respects  as 
are  the  methods  of  analysis  applied  to  them.   The  experimental 
arrangements  can  be  generally  classified  in  two  groups:  (I)  steady- 
state  type  experiments  in  which  flow  streams  of  different  composi- 
tions pass  on  opposite  sides  of  a  barrier  through  which  diffusion 
occurs,  or  (II)  transient  experiments  in  which  diffusion  elimin- 
ates an  initial  concentration  difference  between  two  parts  of  a 
closed  batch  system.   The  methods  of  analysis  may  similarly  be 
divided  as:  (A)  those  based  on  Fick's  first  law  equation,  or  (B) 
those  based  on  his  second  law  equation.   Systems  which  yield 
themselves  to  method  (A)  generally  have  the  advantage  of  much 
simpler  mathematical  descriptions,  but  frequently  are  more  de- 
manding in  experimental  equipment  and  technique.   Type  (I)  ex- 
periments are  illustrative  of  this  point;  they  have  the  advantage 
of  straightforward  analysis  based  on  Fick's  first  law  equation, 
but  require  large  amounts  of  material  and  precise  flow  control. 
Early  measurements  by  method  (I)  were  made  by  Buckingham  (10)  on 
porous  soil;  more  recent  measurements  on  catalysts  and  absorbents 
include   those  of  Scott  and  Dullien  (67) ,  Otani  and  co-workers 
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(56) ,  Evans  and  co-workers  (22) ,  and  Wakao  and  Smith  (77) .   Type 
(II)  experiments  regardless  of  the  method  of  analysis  require 
measurements  of  concentrations  as  functions  of  time,  and  these 
measurements  have  been  their  primary  source  of  difficulty.   The 
development  of  radioactive  tracer  techniques  and  particularly  the 
development  of  semiconductor  radiation  detectors  has  lessened  the 
difficulty  of  making  internal  measurements  of  concentration  without 
disturbing  the  system. 

The  lack  of  attention  to  the  determination  of  diffusion  coef- 
ficients in  the  transition  pressure  range  has  been  notable  in  the 
past.   Only  recently  have  attempts  been  made  to  measure  gaseous 
mutual  and  self-diffusion  coefficients  under  conditions  of  mixed 
Knudsen  and  bulk  diffusive  flow.   The  apparatus  used  in  this  experi- 
ment was  designed  to  measure  diffusive  flow  through  a  capillary  tube 
under  these  conditions  by  the  use  of  radioactive  tracers  and  a  semi- 
conductor radiation  detector.   While  other  investigators  have  uti- 
lized radioactive  tracer  techniques  it  is  believed  that  Mistler's 
(5  3)  iii  the  only  other  work  using  a  semiconductor  detector.   The 
following  survey  of  earlier  techniques  serves  to  illustrate  the  ad- 
vantages of  this  method. 

2.2   Methods  of  Measuring  Diffusion  Coefficients 

Most  of  the  measurements  of  the  gaseous  mutual-diffusion  coef- 
ficient have  been  based  on  a  method  attributed  to  Loschmidt  (51) . 
In  its  simplest  form,  the  apparatus  consists  of  a  long,  closed, 
vertical  tube  that  can  be  divided  into  two  equal  sections  by  a 
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removable  partition  at  its  center.   The  upper  half  is  filled  with 
the  lighter  gas  and  the  lower  half  with  the  heavier  gas,  both  at 
the  same  temperature  and  pressure.   The  partition  is  removed  and 
the  two  gases  are  allowed  to  diffuse  into  each  other  for  a  measured 
amount  of  time. 

Concentration  measurements  may  be  in  the  form  of  a  detection 
device,  generally  at  one  end  of  the  tube,  or  the  total  amount  of  the 
diffusing  material  in  one  of  the  cylinders  may  be  determined  by 
suitable  means.   Either  method  requires  solution  of  the  partial 
differential  equation  produced  by  application  of  Fick's  second  law 
equation.   (Jost  (44)  and  Crank  (15)  give  rather  complete  summaries 
of  the  mathematical  description  of  the  various  geometrical  diffusion 
assemblies  requiring  analysis  by  method  (B)  with  various  initial 
concentration  distributions) .   In  the  case  of  concentration  measure- 
ments by  use  of  a  detector,  the  Fourier  power  series  solution  of  the 
concentration,  in  terms  of  the  diffusion  coefficient,  tube  length, 
and  time,  may  be  used  to  evaluate  the  diffusion  coefficient.   In 
this  latter  case,  the  data,  consisting  of  the  time  for  diffusion 
and  the  composition  of  the  gas  mixtures  contained  in  both  portions 
of  the  cylinder,  are  put  into  proper  forms  of  the  integrated  Fourier 
power  series  and  the  diffusion  coefficient  is  evaluated.   Tordai 
(74)  has  presented  an  analysis  of  the  Loschmidt  method  showing  the 
possibility  of  accuracies  to  within  0.3%.   Strehlow  (71)  has  re- 
ported values  with  less  than  2%  error  by  use  of  a  modified  Loschmidt 
cell. 

An  interferometric  method  of  measuring  mutual  diffusion  coef- 
ficients based  on  the  Loschmidt  geometry  was  presented  by  Boyd, 


20 


Stein,  Steingrimsson,  and  Rumpel  (7)  in  1951.   In  this  method  the 
change  in  composition  of  the  gas  at  a  fixed  position  in  the  dif- 
fusion cell  is  determined  by  measuring  the  change  in  refractive 
index  with  an  interferometer.   Several  determinations  of  the  bi- 
nary diffusion  coefficient,  D. 2 ,  can  be  made  during  the  course  of 
a  run,  and  values  accurate  to  within  0.5%  have  been  obtained  by 
this  method. 

Boardman  and  Wild  (4)  simulated  the  conditions  for  self- 
diffusion  by  using  gases  with  molecules  of  equal  mass  in  a 
Loschmidt  cell.   Using  mixtures  of  nitrous  oxide  and  carbon  di- 
oxide, D-^  for  either  gas  was  found  to  be  0.0107  cm2/sec  at  15°C 
and  one  atmosphere  pressure.   By  the  method  of  Kelvin  triads,  (45) 
they  also  determined  D-i-,  for  hydrogen,  nitrogen,  and  again  for 
carbon  dioxide.   In  this  method,  D12  is  determined  for  each  of 
the  three  pairs  of  gases  found  in  a  ternary  system  and  o12  is 
calculated.   Since  o-.,  =  (1/2)  (oil  +  "22'  '  values  f°r  eacn  of 
the  three  0,2's  allow  solution  for  each  value  of  o^^. 

The  Loschmidt  method  and  related  methods  are  quite  time  con- 
suming, and  it  is  very  difficult  to  estimate  the  effect  of  con- 
centration on  diffusivity.   Coward  (14)  gives  an  indication  of 
the  time  required  for  Loschmidt  type  determinations.   He  indicated 
that  each  run  is  allowed  to  go  on  for  sixty  minutes  and  then  the 
mixtures  are  allowed  to  stand  overnight  before  being  analyzed. 
Boardman  (4)  indicates  minimum  diffusion  time  of  one  hour  and  a 
maximum  of  six  hours. 
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The  method  of  Loschmidt  is  quite  inconvenient  and  difficult 
to  use  in  evaluating  the  diffusion  coefficient  of  liquid  vapors 
into  gases.   Stefan  (70)  devised  a  method  which  makes  such  eval- 
uations relatively  simple  and  convenient  for  reasonably  volatile 
liquids.   A  sample  of  the  liquid,  whose  vapor  diffusion  coeffi- 
cient is  to  be  measured,  is  placed  in  the  bottom  of  a  vertical 
tube.   With  the  tube  maintained  at  constant  temperature,  the  sec- 
ond component,  a  gas,  is  passed  over  the  top  of  the  tube  at  a 
rate  sufficient  to  keep  the  partial  pressure  of  the  vapor  there 
at  a  value  essentially  corresponding  to  the  initial  composition 
of  the  gas,  i.e.,  zero  in  the  case  of  a  vapor-free  gas.   The  dif- 
fusion of  the  vapor  out  of  the  tube  affects  a  slow  lowering  of 
liquid  level.   The  rate  of  evaporation  of  the  liquid  is  simply 
related  to  the  diffusion  coefficient.   As  with  the  Loschmidt 
method,  the  data,  consisting  of  the  time  of  gas  flow  and  the 
amount  of  liquid  evaporated,  are  put  into  proper  integrated  forms 
of  the  diffusion  equations  and  D.,  evaluated. 

In  1957  Walker  and  Westenberg  (78)  introduced  a  significant 
new  flow  method  for  measuring  the  diffusion  coefficient  of  gases 
called  "the  point  source  technique".   This  shortened  significantly 
the  time  period  required  for  an  experiment.   The  method  made  use 
of  the  steady  injection  of  a  trace  gas  from  a  fine  hypodermic 
tube  into  a  slow,  laminar,  uniform  stream  of  a  second  carrier  gas. 
Measurements  were  made  by  thermal  conductivity  cells  on  the  con- 
centration of  the  trace  gas  in  both  the  radial  and  longitudinal 
directions.   The  diffusion  coefficient  could  be  calculated  from 
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either  the  radial  or  longitudinal  data  with  an  error  of  approxi- 
mately 1%.   By  heating  the  carrier  gas  the  measurements  may  be 
extended  to  fairly  high  temperatures.   Also  the  trace  gas  concen- 
trations involved  are  essentially  zero,  so  that  uncertainties  re- 
garding the  concentration  dependence  of  the  diffusion  coefficient 
are  avoided.   The  method  has  been  used  to  study  He-N2 ,  N2-C02, 
He-A,  and  C02~He-N2  (78,  79,  80,  81). 

Ember,  Ferron,  and  Wohl  (20)  have  modified  the  point  source 
technique  to  permit  investigation  at  higher  temperatures  through 
the  use  of  a  thin,  flat  flame  which  burns  at  the  base  of  the 
apparatus.   They  also  used  radioactive  tracer  techniques  enabling 
them  to  measure  self-diffusion  coefficients. 

Bohemen  and  Purnell  (5)  and  Giddings  and  Seager  (27,  28,  29) 
first  successfully  demonstrated  that  diffusion  coefficients  at 
low  pressures  could  be  measured  accurately  by  the  use  of  gas 
chromatographic  techniques.   The  apparatus  used  consisted  of  a 
gas  chromatography  unit  with  an  empty  tube  replacing  the  packed 
column.   Reported  binary  coefficients  agreed  with  published  re- 
sults within  1%.   Giddings   et  al  (30)  and  Evans  and  Kennedy  (21) 
have  extended  this  technique  over  a  range  of  temperatures  and 
pressures.   Knox  and  McLaren  (46)  modified  the  method  to  another 
method  called  "arrested  elution  method"  for  measurement  of  diffu- 
sion coefficients  at  low  pressures.   They  reported  D12  =  0.165 
cm2/sec  for  ethylene  in  nitrogen  at  18°C  and  750  mm  which  was  in 
agreement  with  the  literature  value. 
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Anudsen  (47) ,  Gaede  (26) ,  and  Adzumi  (1)  were  the  first  to 
investigate  the  flow  of  gases  in  the  transition  pressure  range. 
Knudsen's  apparatus  consisted  of  two  flasks  of  approximately  1000 
ml  volume  each,  joined  by  a  capillary  tube  with  a  mercury  cutoff 
valve.   He  found  that  a  minimum  specific  flow  occurred  at  a 
pressure  corresponding  to  approximately  r/n   =  0.3.   Gaede  re- 
peated his  measurements  with  an  arrangement  which  used  a  rectan- 
gular slit  in  place  of  the  capillary  tube.   He  verified  Knudsen's 
findings  but  erroneously  attributed  the  minimum  to  an  absorbed 
layer  of  molecules  which  restricted  flow.   Adzumi  also  verified 
Knudsen's  results  with  similar  equipment  and  used  several  tubes 
of  different  bores  at  both  low  and  intermediate  pressures. 

The  closest  approximation  to  self-diffusion  is  accomplished 
by  use  of  isotopic  molecules.   The  intermolecular  forces  and 
effective  cross  sections  for  interactions  between  isotopic  mole- 
cules are  essentially  the  same  as  for  interactions  between  iden- 
tical molecules  of  the  isotopic  pair.   The  only  noticeable  dif- 
ference between  the  isotopic  diffusion  coefficient  and  the  idea- 
lized self-diffusion  coefficient  results  from  the  change  in  the 
reduced  mass  which  is  M,/2  or  M2/2  for  self-diffusion.   The  first 
isotopic  diffusion  measurements  were  made  by  Harteck  and  Schmidt 
(33)  for  the  interdif fusion  of  the  ortho  and  para  forms  of  hydro- 
gen using  a  Loschmidt  cell.   This  method  was  later  used  by  Groth 
and  Harteck  to  determine  the  coefficients  of  krypton  and  xenon 
(32) .   A  more  accurate  method  of  making  isotopic-dif fusion  mea- 
surements was  developed  in  1943  by  Ney  and  Armistead  (54) .   Since 
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nearly  all  subsequent  measurements  of  self-diffusion  have  utilized 

the  Ney-Armistead  technique,  and  since  it  is  most  similar  to  that 

used  in  this  study,  it  will  be  described  in  some  detail. 

The  Ney-Armistead  apparatus  differs  from  that  of  Loschmidt 

in  that  two  large  chambers  are  used  with  a  straight  tube  or  pipe 

connecting  them.   In  the  original  apparatus  the  connecting  tube 

consisted  of  two  eaual-sized  copper  pipes  with  a  section  of  neo- 

prene  tubing  between  them;  by  pinching  the  rubber  with  a  clamp, 

the  two-chambers  could  be  isolated.   One  chamber  was  filled  with 

a  sample  enriched  in  U2   F,  and  the  other  with  a  sample  enriched 

in  U238F, .   With  both  chambers  at  the  same  pressure  and  tempera- 
6 

ture,  the  diffusion  process  was  started  by  unpinching  the  neoprene. 
Samples  from  one  chamber  were  slowly  bled  through  an  adjustable 
capillary  leak  into  the  ion  source  of  a  mass  spectrometer  in  order 
to  follow  the  concentration  changes.   Diffusion  coefficients  were 
reported  to  within  4%.   This  method  has  two  advantages  over  the 
Loschmidt  technique:  (1)  the  geometry  is  more  compact,  reducing 
the  errors  associated  with  thermal  gradient;  and  (2)  continuous 
concentration  measurements  allow  several  determinations  of  the 
diffusion  coefficient  within  one  experiment. 

The  theory  of  the  Ney-Armistead  method  is  based  on  four 
assumptions:  (1)  that  the  diffusion  coefficient  is  independent 
of  concentration;  (2)  that  the  composition  of  the  gas  in  the 
chambers  is  uniform,  the  concentration  gradient  being  confined  to 
the  connecting  tube;  (3)  that  the  concentration  in  the  connecting 
tube  varies  linearly  with  distance  along  the  tube;  and  (4)  that 
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the  volume  of  the  connecting  tube  is  negligible  compared  with  the 
volume  of  the  chambers.   Assumption  (1)  is  justified  for  self- 
diffusion  but  not  for  mutual  diffusion;  according  to  kinetic  theory 
the  mutual  diffusion  coefficient  is  independent  of  concentration 
only  to  a  first  approximation.   A  discussion  of  the  validity  of 
the  other  assumptions  has  been  presented  by  Paul  (58) .   The  sec- 
ond assumption  is  justified  providing  a  small  end  correction  is 
made.   Actually,  the  concentration  gradient  extends  slightly  be- 
yond the  ends  of  the  tube,  and  the  correction  is  made  by  replacing 
the  measured  length  of  the  tube  by  a  slightly  greater  effective 
length.   Analysis  of  the  analogous  electrical  case  by  Maxwell 
(52)  and  Rayleigh  (63)  shows  that  the  length  of  the  tube  should 
be  replaced  by  an  effective  length  1.642r  greater.   Assumption 
(3)  is  found  to  be  valid  after  the  decay  of  initial  transients. 
This  decay  requires  only  a  few  seconds,  a  negligible  time  in  the 
original  experiments.   Finally,  a  small  correction  proportional 
to  the  ratio  of  the  volume  of  the  capillary  tube  to  that  of  the 
chambers  may  be  made  in  order  to  account  for  the  effect  of  the 
tube  volume.   Based  on  the  above  assumptions,  Ney  and  Armistead 
applied  Fick ' s  first  law  to  obtain  an  expression  for  the  concen- 
tration of  an  isotope  in  the  chamber  being  analyzed  as  a  function 
of  time.   Knowledge  of  the  physical  parameters  allowed  a  solution 
for  »    diffusion  coefficient. 

Hutchinson  (39)  determined  the  self-diffusion  coefficient 
of  argon  over  a  range  of  temperatures  by  observing  the  diffusion 
of  A   with  a  modified  Ney-Armistead  apparatus.   He  showed  that 
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the  self-diffusion  coefficient  of  a  gas  can  be  obtained  from  the 
coefficient  of  diffusion  of  one  isotope  into  another  by  the 
following  relationship: 

Du  -  [2m2/(m1+m2)]1/2D12  (25) 

where  m,  and  m,  refer  to  the  molecular  masses  of  the  isotopes. 
The  concentrations  of  A   in  both  chambers  were  measured  by  the 
ion  currents  drawn  to  insulated  electrodes  inserted  into  each 
chamber.   Winn  investigated  the  N14N14— N1  N15  system  (84)  and 
the  temperature  dependence  of  the  self-diffusion  coefficients 
of  A,  Ne,  N2,  02,  C02,  and  CH4  (85)  using  the  Ney-Armistead 
apparatus.   Winn  and  Ney  (83)  used  the  same  method  for  investi- 
gation of  CH.  self-diffusion  for  pressures  ranging  from  515  mm 
to  630  mm  and  temperatures  ranging  from  19.2°C  to  24.7°C. 

The  radioactive  tracer,  C1402  has  been  frequently  used  to 
measure  diffusion  coefficients.   Admur  (2)  measured  diffusion 
coefficients  within  2%  for  C02"C02,  C02-N20,Xc-Xe ,  and  ?-Xe  sys- 
tems using  a  Loschmidt  device.   The  ionization  currents  from 
electrodes  in  the  cells  were  used  to  measure  concentration  changes. 
Drickamer  (40,73)  utilized  a  scintillation  detector  in  a  modified 
Loschmidt  cell  to  measure  diffusion  coefficients  of  C02-C02  and 
C0,-CH.  systems  over  a  pressure  range  from  0 . 5  to  several  atmo- 
spheres.  Diffusion  coefficients  for  C02  were  determined  by  O'Hern 
and  Martin  (55)  from  0°  to  100°C  and  at  pressures  to  200  atm. 
They  used  a  modified  Ney-Armistead  apparatus.   Winter  (86,87)  used 
a  Nev-Armistead  arrangement  to  measure  diffusion  coefficients  of 
C02,  02,  and  N2 .   Visner  (75)  has  made  the  only  study  of  gas 
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diffusion  in  the  transition  pressure  range  using  a  diffusion  tube 
with  well  defined  dimensions,  i.e.,  an  empty  capillary  tube  in- 
stead of  a  chamber  packed  with  porous  material.   He  used  a  modi- 
fied Mey-Armistead  cell  with  end-window  Geiger  counters  installed 
in  both  flasks  to  follow  concentration  changes.   The  self-diffusion 
coefficient  of  xenon  was  measured  for  capillaries  of  0.0025  cm  and 
0.025  cm  radius  in  tubes  of  three  lengths  for  pressures  corre- 
sponding to  radius  to  mean  free  path  ratios  from  0.001  to  65. 
The  trace  gas  was  injected  into  one  flask  after  filling  the  flasks 
to  the  desired  pressure  with  normal  gas  by  breaking  a  tiny  glass 
bulb  of  gas  at  the  correct  pressure  to  prevent  the  creation  of  a 
pressure  gradient  across  the  capillary  tube. 

An  apparatus  called  a  diffusion  bridge  was  used  by  Bendt  (3) 

3    4 
to  measure  the  diffusion  coefficient  for  He  -He   and  H2-D2  at  one 

atmosphere  pressure  and  temperatures  from  1.74°K  to  296°K.   The 
measurement  was  made  while  steady-state  flow  existed  through  four 
capillaries  and  through  the  diffusion  tube  which  joined  the  supply 
and  exhaust  capillaries  in  the  temperature  bath.   The  pressure 
gradient  in  the  diffusion  tube  was  made  as  small  as  possible  by 
adjusting  the  flow  rates  of  the  gases  through  the  capillary  feed 
lines.   The  exhaust  capillaries  were  alternately  connected  to  a 
mass  spectrometer  for  gas  analysis.   Small  fluctuations  in  the 
flow  rates  and  in  the  sensitivity  of  the  mass  spectrometer  limited 
the  accuracy  of  the  measurements  to  an  error  of  2-6%.   The  cor- 
rection factors  used  were  the  Rayleigh  tube  length  correction  and 
the  mass  correction  formula  (Eq.  25)  for  isotopic  self-diffusion 
measurements . 
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14 
In  contrast  to  the  numerous  reports  of  use  of  C  Oj  as  a 

tracer  in  diffusion  measurements,  no  mention  of  a  similar  use  of 

C14H.  was  found  in  the  literature.   Winn  (85)  and  Winn  and  Ney 

(83)  used  C13H  -C12H.  for  isotopic  diffusion  measurements  of  CH .  , 
v  4      4 

but  the  stability  of  C13  precluded  the  use  of  tracer  techniques 

and  they  relied  upon  mass  soectrography  for  concentration  measure- 

14 
merits .   It  appears  quite  possible  that  the  use  of  C  H4  as  a 

tracer  in  determination  of  self-diffusion  coefficients  for  CH4  is 

a  unique  feature  of  this  study. 

The  above  survey  has  included  a  discussion  of  all  common 

techniques  used  to  measure  diffusion  coefficients.   The  Ney- 

Armistead  method  has  been  emphasized  because  it  allows  a  simple 

and  accurate  means  of  studying  gas  diffusion  in  the  transition 

pressure  range.   Additional  surveys  of  recent  experimental  and 

theoretical  work  on  the  diffusion  of  gases  have  been  made  by 

Liley  (50)  and  Johnson  (41,42,43).   A  discussion  of  information 

on  diffusivity  needed  by  the  chemical  and  process  industries  has 

been  presented  by  Friend  and  Adler  (25) . 

2.3   Analysis 

The  diffusion  cell  used  in  this  study  was  of  the  Ney- 
Armistead  type  and  is  depicted  in  Fig.  2.   To  make  a  measurement 
of  the  effective  diffusion  coefficient,  both  flasks  were  evacuated 
and  then  filled  with  gas  at  the  same  temperature  and  pressure. 
The  gas  in  flask  I  included  a  trace  of  radioactive  gas.   (See 
Appendix  C  for  a  detailed  description  of  experimental  procedure) . 
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Following  the  opening  of  the  capillary  tube  valve  a  one  dimen- 
sional diffusion  process  was  established  in  which  the  molecules 
of  radioactive  gas  from  flask  I  intermixed  with  the  gas  from  flask 
II  until  the  system  was  of  constant  composition  throughout  its 
volume.   The  amount  of  material  transferred  by  self-diffusion  was 
then  directly  proportional  to  the  change  in  radioactivity  in  flask 
I.   This  change  in  activity  as  a  function  of  time  was  measured  by 
use  of  a  semiconductor  radiation  detector. 

Description  of  the  diffusion  process  in  the  tube  by  Fick's 
second  law  caused  only  Ney  and  Armistead's  first  and  second  assump- 
tion to  be  required.   Based  upon  these  assumptions  and  the  use  of 
Fick's  second  law  the  resulting  partial  differential  equation  for 
diffusion  through  the  tube  was 

|§-  D0  (26) 

where  C  =  the  concentration  of  radioactive  gas  at  time  t  at  posi- 
tion z,  D  =  the  diffusion  coefficient,  and  z  =  the  distance  along 
the  tube  length. 

The  boundary  conditions  for  Eq.  (26)  were: 

t  =  0,  z  >  0,  C  =  Cj  (27) 


t  >  0,  z  =  0,  C|z=f)  =  C0  -  V"1  f  (-D„R2)ff|    dt 

U  i  2=0 


(28) 
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9C 

3z| 


t    >    0,    z   =   L,    C|z=L   =   C,    +    V-^JVd.R2)- 

f    TtR2Cdz    +    C.irR2L)     =    C,    +    vl1     f     (-DirR2)|£  dt 

Jn  X  1  2      Jo  3Z|z=L 


(29) 


where  V.  =  volume  of  flask  I  (tagged) 

V?   =  volume  of  flask  II  (untagged) 

L  =  effective  tube  length 

R  =  radius  of  tube 

C   =  initial  concentration  of  radioactive  gas  in  flask  I 
o 

C,  =  initial  concentration  of  radioactive  gas  in  flask  II 

and  tube 
C  =  concentration  of  radioactive  gas  as  a  function  of  t 

and  z . 
Eciuation  (27)  states  that  at  time  zero  the  concentration  of  tagged 
material  in  the  tube  and  flask  II  is  C. .   Equation  (28)  shows 
that  the  amount  of  radioactive  gas  in  flask  I  is  the  original 
amount  less  that  leaving  the  flask  by  diffusion  into  the  capillary 
tube.   The  final  boundary  condition,  Eq.  (29)  states  that  the 
amount  of  tracer  material  in  flask  II  is  the  original  amount  plus 
that  entering  from  the  tube;  the  amount  entering  from  the  tube 
may  be  evaluated  directly  as  in  Eq.  29  or  as  the  total  of  that 
leaving  the  flask  I  and  that  originally  in  the  tube  less  the  amount 
remaining  in  the  tube.   This  latter  form  is  shown  in  Eq .  (29) . 
The  complete  solution  to  Eq.  (26)  is 
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C(x,f)    =    (8+aX,  [B+l])  (a+B+c.S)~      +    (X-l)     [     [exp(-X    t)  ] 
1  L        n=l 


[X    cosX   +BsinX    1  [acos(X   e)    -   »    sin(»   {)] 
n  n  n  n'  n  n 


{X2[(2x    )-1(a+B+aB-X2)cosX„   -   sinX 

n  n  n  n  n 


(l+(a+B)2)}    X         .  <3n> 


where    X      are   the   non   zero   roots   of 
n 

tan    X      =    Xn(a+B)  (X^-aB)"1         .  (31) 

A  detailed  solution  of  Eq.  (31)  by  the  method  of  Laplace  trans- 
forms is  given  in  Appendix  A  as  well  as  a  description  of  the 
dimensionless  variables  employed. 

The  time  response  of  the  system  used  is  shown  in  Figs.  (3) 
and  (4)  where  the  concentration  of  tracer  in  flask  I  is  shown  as 
a  function  of  time  for  selected  pressures.   Figure  (3)  is  for 
self-diffusion  in  a  C02  system  while  Fig.  (4)  is  for  a  CH4  system. 
It  should  be  noted  that  extremely  accurate  short  time  measurements 
are  not  required  and  that  the  rate  of  response  may  be  easily 
varied  by  changing  the  tube  length. 

Equation  (30)  relates  the  concentration  of  the  tagged  gas  to 
time  at  a  position  z  inside  the  cell  in  terms  of  the  diffusion 
coefficient,  D,  and  the  cell  dimensions.   An  application  of  least 
squares  analysis  (33)  to  the  concentration  versus  time  data  using 
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Eq.  (30)  and  the  appropriate  system  parameters  yields  an  evalua- 
tion of  the  diffusion  coefficient.   A  detailed  discussion  of  the 
procedure  for  applying  least  squares  analysis  to  the  evaluation 
of  D  is  given  in  Appendix  C  as  well  as  an  estimate  of  the  error 
involved.   In  this  experiment  the  tube  length  was  determined  by 
measuring  and  applying  the  Rayleigh  end-correction  used  by  Ney 
and  Armistead.   The  other  system  parameters  were  measured  as  out- 
lined in  section  2.5. 

Based  upon  the  assumption  of  continuous  mixing  at  the  pres- 
sure of  the  measurement,  which  has  been  shown  to  be  valid  by 
Visner  (74)  ,  the  radioactivity  measured  by  the  surface-barrier 
detector  was  linearly  proportional  to  the  concentration  of  tracer 
in  flask  I.   Because  of  the  rapid  response,  the  thin  window,  and 
the  high  counting  efficiency  of  the  detector  it  was  not  necessary 
to  make  dead  time  or  other  corrections  on  the  activity  data.   The 
properties  of  semiconductor  detectors  are  discussed  in  section 
2.4. 

2.4   Semiconductor  Radiation  Detectors 

As  a  result  of  rapid  technological  advances  in  design  and 
production  technioues  semiconductor  radiation  detectors  have  at- 
tained wide  acceptance  and  usage  in  the  last  few  years  (8,16,17, 
37,38,61).   The  demand  for  their  development  is  readily  explained 
by  a  brief  review  of  their  advantages.   A  semiconductor  radiation 
detector  approximates  an  ideal  ionization  chamber.   This  is  ac- 
complished by  replacing  the  conventional  gas  of  the  ionization 
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chamber  by  solid  semiconductor  materials  of  high  stopping  power. 
The  linear  response  of  the  signal  pulse  size  (1  electron  pair/3.5 
ev)  to  the  energy  dissipated  in  the  detector  by  an  ionizing  parti- 
cle make  semiconductor  detectors  particularly  suitable  for  measure- 
ments of  energy  spectra.   The  high  stopping  power  greatly  reduces 
the  probability  of  a  particle  traversing  the  active  volume  without 
complete  dissipation  of  its  energy.   Semiconductor  detectors  have 
the  further  advantages  of  small  size,  low  weight,  thin  window 
thickness,  low  operating  voltage,  precisely  controlled  sensitive 
depth  and  geometry,  rapid  time  response,  stability,  and  durability. 

Semiconductor  detectors  utilize  the  properties  of  a  junction 
of  an  n-type,  donor-rich  material  with  a  p-type ,  acceptor-rich 
material  (69,72).   The  application  of  a  reverse  electrical  bias 
to  the  junction  causes  a  migration  of  electrons  in  the  n-type 
material  toward  the  p-type  material  and  holes  in  the  p-type  toward 
the  n-type.   The  results  of  these  movements  are  the  establishment 
of  a  space  charge  region  caused  by  an  unbalanced  concentration  of 
filled  electron  acceptor  sites  and  empty  donor  sites  on  opposite 
sides  of  the  junction  and  the  formation  of  a  carrier  free  region 
of  high  resistance.   This  carrier  free  region  is  termed  the  de- 
pletion zone.   Its  depth  is  proportional  to  the  square  root  of  the 
product  of  the  bias  voltage  and  the  dielectric  resistivity.   This 
depth  generally  ranges  from  2  to  3  mm  depending  on  the  type  of 
detector. 

Goulding  (31)  presents  a  discussion  of  the  applications  and 
properties  of  the  different  types  of  semiconductor  detectors. 
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The  three  commonly  used  types  are:  (1)  the  silicon-diffused 
junction  type;  (2)  the  lithium-drifted  silicon  or  germanium  type; 
and  (3)  the  silicon  surface-barrier  type.   The  methods  of  prepar- 
ation of  these  types  in  the  same  order  are:  (1)  diffusion  of  n- 
type  Dhosphorous  into  p-type  (boron  doped)  silicon;  (2)  injection 
of  n-type  lithium  ions  deep  into  p-type  silicon  or  germanium;  and 
(3)  evaporation  of  a  thin  layer  of  p-type  gold  on  the  surface  of 
n-type  silicon. 

A  silicon  surface-barrier  detector  was  used  in  this  study. 
If  an  ionizing  particle  possesses  sufficient  energy  to  penetrate 
the  thin  gold  film  and  a  portion  of  the  depletion  region  and 
impinges  the  surface  at  an  angle  which  allows  this,  it  will  create 
electron  pairs  as  it  dissipates  its  energy  in  the  sensitive  volume. 
The  electron  pairs  are  separated  by  the  space  charge  and  attain  a 
drift  velocity  as  they  migrate  toward  the  electrical  poles  of  the 
space  charge.   This  motion  induces  a  charge  in  the  bias  circuitry. 
By  use  of  a  charge  sensitive  preamplifier  with  negative  capacitive 
feedback  this  charge  may  be  converted  to  a  voltage  pulse  which  is 
essentially  independent  of  detector  capacitance. 

The  accuracy  of  measurements  made  with  a  surface-barrier 
detector  is  directly  related  to  the  ratio  of  the  rate  of  true  sig- 
nal propagation  to  that  of  noise  or  spurious  signals.   Spurious 
signals  may  be  generated  within  the  amplification  system  or  with- 
in the  detector  itself.   Noise  generation  within  the  amplification 
system  is  best  avoided  by  careful  selection  of  the  preamplifier 
and  amplifiers.   Discussions  of  these  considerations  are  presented 
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bv  Fairstein  (24)  and  Radeka  (62) .   The  primary  sources  of  detec- 
tor noise  are  microplasma  breakdown,  edge  leakage  currents,  and 
charge  carrier  generation  and  recombination  in  the  depleted  volume. 
Detector  noise  may  be  minimized  by  selection  of  a  high  quality  de- 
tector with  a  small,  thin  depletion  region  and  operation  at  mini- 
mum temperatures.   The  detector  should  also  have  low  electrical 
capacitance,  but  the  capacitance  is  proportional  to  the  ratio  of 
the  area  of  the  depletion  region  to  its  thickness.   Thus  a  com- 
promise must  be  made  in  selecting  the  depletion  thickness.   The 
choice  depends  largely  upon  the  intended  use. 

The  use  of  a  surface-barrier  detector  requires  that  several 
other  requirements  be  met.   All  signal  leads  and  the  detector 
must  be  shielded  to  prevent  the  generation  of  noise  resulting 
from  electromagnetic  radiations.   The  detector  is  also  sensitive 
to  visible  light  and  must  be  in  a  thoroughly  darkened  enclosure. 
The  signal  lead  from  the  detector  to  preamplifier  should  be  of 
minimum  length  to  reduce  input  capacitance.   All  electronic 
equipment  should  be  connected  to  common  grounds  to  avoid  ground 
loop  problems . 

The  detector  bias  voltage  must  be  stable  and  well  regulated. 
This  voltage  should  be  chosen  so  as  to  give  the  best  signal  to 
noise  ratio.   The  speed  and  efficiency  of  charge  collection  with- 
in the  detector  increase  monotonically  with  increasing  voltage. 
The  noise  level  also  increases  with  detector  bias  following  an 
initial  decrease  until  a  breakdown  of  the  detector  occurs.   The 
noise  level  may  be  monitored  as  bias  voltage  is  increased  by  use 
of  an  RMS  voltmeter,  an  oscilloscope,  or  a  multichannel  analyzer. 
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Every  precaution  should  be  taken  to  avoid  damaging  the  de- 
tector.  The  fragile  gold  surface-barrier  must  not  be  damaged  by 
the  application  of  sudden  or  combined  stresses.   The  detector 
must  be  protected  from  chemical  contamination,  continued  exposure 
to  high  radiation  fluxes,  and  excess  bias  voltage. 

2.5.0  Experimental  Facilities 

A  diffusion  cell  apparatus  and  an  electronic  counting  system 
comprised  the  experimental  facilities  used  in  making  all  measure- 
ments.  With  a  radioactive  component  diffusing  from  one  flask  to 
a  second  flask  joined  by  a  capillary  tube,  the  radioactivity  in 
the  charged  flask  provided  a  continuous  measurement  of  the  diffu- 
sion occurring  through  the  opening  in  the  capillary  tube.   The 
diffusion  cell  employed  has  similarities  to  those  of  earlier  in- 
vestigators, especially  to  those  of  Hutchinson  (39),  O'Hern  and 
Martin  (55)  ,  Winter  (87)  ,  Ney  and  Armistead  (54)  ,  and  Timmerhaus 
and  Drickamer  (73).   Figures  2,  5 ,  6  and  7  illustrate  the  diffusion 
cell  and  auxiliary  apparatus.   The  electronic  counting  system, 
diagrammed  in  Fig.  8  consisted  of  a  semiconductor  radiation  detec- 
tor and  a  scaling  circuit. 

2.5.1   Diffusion  Cell  and  Auxiliary  Apparatus 

A  top  view  of  the  Pyrex  glass  diffusion  cell  mounted  in  the 
unfiLled  water  bath  is  shown  in  Fig.  6.   Referring  to  Fig.  5, 
flask  I  represents  the  side  of  the  system  charged  with  radioactive 
carbon  dioxide  and  containing  the  surface-barrier  semiconductor 
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Figure  7. 


Surface-barrier  radiation  detector  before  insertion 
into  operating  position. 
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Figure  8.   Block  diagram  of  complete  counting  system. 
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radiation  detector.   The  detector  was  installed  in  flask  I,  sus- 
pended from  two  No.  22  tungsten  wire  leads  passing  through  a  re- 
movable S  34/45  ground  glass  stopper.   The  capillary  tube  and  a 
pressure  equalization  tube  "11"  provided  vacuum  tight  connections 
between  flask  I  and  flask  II.   The  length  and  the  internal  dia- 
meter of  the  precision  bore  capillary  tube  were  found  to  be 
10.00   cm  and  0.1026  cm  respectively.   The  diameter  was  determined 
by  measuring  the  length  and  mass  of  a  thread  of  mercury  inside  the 
tube.   High  vacuum  stopcocks  were  used  to  control  gas  flow  through 
the  six  orifices  in  flask  I  and  II,  through  the  pressure  equali- 
zation tube,  and  to  the  thermocouple  vacuum  gauge  tubes,  "6"  and 
"18  '  (See  Appendix  B  for  detailed  experimental  procedure)  .   Flasks 
I  and  II  had  volumes  of  579  and  512  ml,  respectively.   The  capaci- 
ties of  these  flasks,  up  to  the  bottom  of  the  ground  glass  joint 
(see  Fig.  2) ,  were  determined  by  filling  with  a  measured  quantity 
of  water. 

A  stiff  molybdenum  wire  passing  through  a  spherical  ground 
glass  joint  provided  a  lever  arm  for  operation  of  the  capillary 
tube  valve  "10".   A  small  conical  piece  of  soft  rubber  was  attached 
to  the  end  of  the  wire  so  that  it  entered  and  sealed  the  end  of 
the  capillary  tube  when  the  valve  was  in  the  closed  position. 
This  capillary  tube  valve,  which  differed  from  those  of  many 
earlier  experimental  diffusion  arrangements,  was  designed  to  pre- 
vent any  discontinuities  or  irregularities  along  the  diffusion 
path. 
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Radioactive  and  nonradioactive  gas  flowed  to  flask  I  through 
stopcocks  "4"  and  "9"  respectively.   The  nonradioactive  gas 
(greater  than  99.5%  pure)  was  purified  in  a  tube  of  Drierite,  silica 
gel,  and  magnesium  perchlorate.   Volk  Radiochemical  Company  sup- 
plied glass  ampoules  of  both  C1402  and  C14H4 .   The  specific  acti- 
vity of  the  C1402  was  1  millicurie/millimole  (0.01%  vol  C   02)  and 
the  ampoule  contained  0.10  mc  of  C   02>   The  total  activity  of 
the  ampoule  of  C14II4  was  0.5  mc  of  C  H4  and  the  specific  activity 
was  5  mc/mm  (0.05%  vol  C14H4> .   After  the  diffusion  cell  was 
purged  with  nonradioactive  gas,  trace  gas  was  throttled  into 
flask  I  through  the  Hoke  needle  valve  "3". 

A  two  station  Cenco  vacuum  gauge  (Type  GMA-140)  was  utilized 
in  measuring  pressures  in  the  flasks.   By  closing  "11",  the  pres- 
sure in  flask  I  was  measured  by  gauge  tube  "6"  with  valve  "5" 
open;  the  pressure  in  flask  II  was  measured  by  gauge  tube  "18" 
with  "11"  closed  and  "13"  open. 

An  upright  McLeod  vacuum  gauge,  "14",  (range  0.05  to  500 
microns)  connected  to  flask  II  was  used  in  calibrating  the  thermo- 
couple gauge.   For  methane  quadratic  relationships  were  found  for 
the  calibrations  between  gauge  readings  and  absolute  pressure  in 
the  range  of  interest  (0-500  microns) .   The  vacuum  gauge  meter 
could  be  read  to  approximately  +1%  between  0  and  100  microns  and 
to  +5%  at  higher  pressures. 

A  12  gallon  tank  filled  with  water  served  as  a  constant 
temperature  (20°C)  bath  for  the  diffusion  cell  (95%  immersed) . 
Constant  temperature  regulation  to  +0.1°C.  was  provided  by  a 
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thermoregulator  and  relay  which  controlled  a  low  wattage  heating 
element.   Tap  water  (16.5°C)  running  through  a  coiled  copper 
tubing  served  as  the  heat  sink. 

The  diffusion  cell,  together  with  all  auxiliary  equipment 
containing  radioactive  gas,  was  housed  in  a  high  velocity  hood  to 
guard  against  exposure  to  radiation.   A  liquid  nitrogen  trap  in- 
stalled between  stopcock  "13"  and  the  two-stage,  oil  sealed 
mechanical  vacuum  pump  prevented  the  escape  of  radioactive  gas 
through  the  pump  exhaust.   The  cold  trap  was  periodically  removed 
to  allow  accumulated  solids  to  sublime  and  be  expelled  slowly 
through  the  hood  exhaust. 

2.5.2   Electronic  Counting  System 

The  ability  of  an  electronic  counting  system  to  detect  low 
energy  beta  particles,  such  as  those  emitted  by  carbon-14 

(E    -  156  kev  (40))  has  been  shown  to  be  dependent  upon  the 

max 

signal  to  noise  ratio  of  the  amplified  detector  output.   If  noise 
amplitude  is  high  compared  to  signal  amplitude,  it  is  usually  not 
possible  to  separate  the  spurious  counts  from  those  that  are  real. 
If,  on  the  other  hand,  the  signal-to-noise  ratio  is  large  it  is 
possible  to  electronically  discriminate  against  the  noise  pulses 

and  count  only  the  pulses  produced  by  radiation  incident  upon  the 

14 
detector.   For  the  case  of  polyenergetic  C   beta  particle  detec- 
tion, part  of  the  real  signal  is  always  discriminated  out  with 
the  noise.   Thus,  the  detection  system  used  in  this  experiment 
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was  designed  to  develop  a  minimum  amount  of  noise.   It  is  most 
important  that  the  production  of  noise  is  prevented  at  that  point 
in  the  counting  system  where  the  real  signal  is  the  smallest,  i.e., 
at  the  detector  where  the  signal  is  created. 

Because  the  surface-barrier  radiation  detector  was  extremely 
sensitive  to  electromagnetic  radiation,  it  was  necessary  to  con- 
struct an  electrical  and  a  light  shield  surrounding  the  glass 
flask  housing  the  detector.   A  fine-mesh  brass  screen  jacket  built 
around  flask  I  (see  Fig.  5)  and  grounded  to  the  preamplifier 
chassis  was  effective  in  shielding  against  stray  electrical  sig- 
nals.  The  complete  outer  surface  of  flask  I,  except  for  small 
parts  of  stopcocks  "4",  "5",  and  "9",  was  painted  with  black 
Glyptal  paint.   Since  the  capillary  tube  served  as  a  light  pipe 
between  the  two  flasks,  part  of  flask  II  was  also  painted.   De- 
tector surface  leakage  was  avoided  by  handling  the  detector  with 
extreme  care.   Thermal  noise  in  the  detector  was  minimized  by 
maintaining  the  detector  at  a  constant  temperature. 

A  five  inch  long  microdot  shielded  cable  connected  the  de- 
tector leads  to  an  ORTEC  low  noise  preamplifier  (Model  101) . 
This  cable  was  made  as  short  as  possible  to  reduce  input  capaci- 
tance since  the  signal  voltage  was  inversely  proportional  to  the 
sum  of  the  effective  detector  capacitance  and  all  other  input 
capacitances.   The  preamplifier  utilized  a  charge-sensitive  input 
circuit  and  a  voltage  amplifier  provided  with  a  negative  feedback 
path  allowing  it  to  simulate  a  charge-to-voltage  transducer. 


48 


Signals  emerging  from  the  preamplifier  had  an  amplitude  of  sever- 
al millivolts. 

The  preamplifier  was  connected  to  the  ORTEC  low  noise  ampli- 
fier (Model  201)  through  a  shielded  cable.   This  amplifier  chassis 
contained  the  power  supplies,  the  60  cycle  test  pulse  generator, 
the  discriminator  bias  circuit,  and  the  main  and  post  amplifiers. 
A  Tektronix  oscilloscope  (Type  515A)  was  used  to  monitor  the  ampli- 
fied signal  and  thereby  select  proper  adjustments  for  the  ampli- 
fiers.  The  positive  signal  output  of  the  post  amplifier  was  in- 
verted and  attenuated  by  a  EG  S  G  IT  100  inverting  transformer  to 
make  it  compatible  with  the  TMC  Gammascope  II  (Model  102)  multi- 
channel analyzer  (scaler  mode) . 

A  square  wave  from  a  Krohn-Hite  low  frequency  generator 
(Model  400-C)  was  used  as  the  external  time  base  for  the  analyzer. 
The  wave  frequency  was  variable  between  0.009  and  1100  cycles  per 
second.   The  paper  tape  containing  the  experimental  data  was 
printed  with  a  Hewlett-Packard  digital  recorder  (Model  561-B) . 
Serious  ground  loop  and  line  transient  problems  were  alleviated 
by  connecting  all  electronic  eouipment,  excluding  the  scope,  to 
a  well  regulated  Beckman  power  supply. 

2.6.0   Reduction  of  Data 

2.6.1   Preparation  of  Data  for  Analysis 

The  method  of  analysis  of  data  in  this  study  was  basically 
the  same  as  that  presented  by  Mistier  (53)  .   It  consists  of  an 
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iterative  least  squares  fitting  of  the  data  to  determine  the 
effective  diffusion  coefficient.   A  description  of  this  analysis 
and  the  computer  program  used  to  perform  it  is  given  in  Appendix 

C. 

The  only  steps  required  to  ready  the  data  for  computer  treat- 
ment were  determination  of  the  average  initial  and  final  count 
rates  and  the  deletion  of  data  points  which  were  obviously  in 
error.   The  average  initial  count  rate  was  found  by  simply  aver- 
aging the  data  points  taken  before  opening  the  capillary  tube 
valve,  while  the  final  average  count  rate  was  the  average  of  data 
points  taken  after  the  system  reached  equilibrium.   The  deletion 
of  some  data  points  within  a  run  was  sometimes  necessitated  be- 
cause the  values  for  these  points  were  obviously  much  too  large. 
The  gross  errors  of  these  values  were  caused  by  vibration  of  the 
system  and  detector  upon  opening  the  capillary  valve  or  by  the 
electronic  pickup  of  electromagnetic  noise.   For  the  data  of  a 
run  to  be  considered  reliable  the  number  of  such  points  encountered 
was  severely  restricted.   A  sample  set  of  data  taken  at  a  pressure 
of  203  u  is  shown  in  Table  I.  An  application  of  least  squares 

analysis  and  the  mass  correction  for  isotopic  self-diffusion  to 

2 
this  data  indicated  a  value  of  Deu  =  614.3  +14.4  cm  /sec.  It  is 

noted  that  the  diffusion  coefficient  actually  measured  was  the 

mutual  diffusion  coefficient  for  the  two  isotopes.   Throughout 

the  analysis,  however,  this  was  assumed  to  be  equal  to  the  self- 

dif fusion  coefficient  of  the  natural  gas.   The  mass  correction 

factor  was  applied  as  a  final  step. 
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Table  I 

Sample  set  of  count-rates  versus 

time  data 

taken  from  the  analyzer  output. 

Timea 

Time 

(sec) 

Counts/Channel        (sec) 

Counts/Channel 

0 

1693              754 

1106 

29 

1654               774 

1069 

49 

1667               793 

1097 

69 

1652               813 

1061 

88 

1561               832 

1096 

108 

1516               852 

1054 

127 

1490               871 

1112 

147 

1483               891 

1057 

166 

1562               911 

1077 

186 

1547               930 

1148 

206 

1505               950 

1036 

225 

1398               969 

1010 

245 

1392               989 

1057 

246 

1378              1028 

1059 

284 

1405              1048 

1042 

304 

1246              1067 

1032 

323 

1481              1087 

1048 

343 

1361              1106 

1021 

362 

1316              1126 

1125 

382 

1272              1148 

951 

401 

1289              1165 

1048 

421 

1247              1185 

9  70 

441 

1288              1204 

1006 

460 

1297              1224 

1004 

480 

1210              1244 

1008 

499 

1266              1263 

1000 

519 

1202              1283 

955 

539 

1156              1302 

1000 

558  . 

1220              1322 

970 

578 

1181              1341 

925 

597 

1147              1361 

993 

617 

1209              1381 

973 

636 

1177              1400 

973 

656 

1167              1420 

976 

676 

1111              1439 

1011 

695 

1209              1459 

949 

715 

1166              1479 

947 

734 

1139              1498 

956 
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Table  I  (Continued) 


Time  Time 

(sec) Counts/Channel [sec) Counts/Channel 

1513  918  1635  955 

1537  1003  1655  974 

1557  976  1674  957 

1576  887  1694  941 

1596  920  1714  896 

1616  913  1733  967 


P  =  203  u 

Initial  pressure  difference  =  -2  y 

D    before  pressure-induced  flow  correction  =  614.3 

ell 
Deviation  arising  only  from  least  squares  analysis  =  12.1 

Combined  deviation  caused  by  least  squares  and  uncertainty  of 
pressure  measurement  =  14.4 

D    after  correction  for  pressure-induced  flow  =  617.4 
ell 


aThe  measured  time  increment  was  19.583  sec;  all  the  times 
above  have  been  rounded  to  the  nearest  second. 
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2.6.2   Standard  Deviation  Analysis 

For  brevity,  the  effective  mutual  isotopic  diffusion  coef- 
ficient, Dei2,  is  referred  to  merely  as  "D"  in  this  section.   The 
uncertainty  associated  with  the  determination  of  the  diffusion 
coefficient  was  the  result  of  scatter  in  the  measured  count  rates, 
a  negligible  error  in  time  measurements  and  uncertainty  in  measure- 
ments of  pressure.   Assuming  a  normal  distribution  for  errors  in 
the  measurement  of  count  rates  an  expression  for  this  uncertainty 
in  D  may  be  found. 

Given  the  relationship 

X  =  f(t,D)  (32) 

where  D  is  a  dependent  variable  of  the  measured  quantities  X  and 
t,  the  error  of  D  due  to  errors  in  measurement  is  given  by  the 
theory  of  propagation  of  errors.   Letting  D  equal  diffusion  coef- 
ficient, t  equal  time,  X  equal  normalized  count  rate,  and  variances 

2     2         2 

of  the  quantities  X,  D,  and  t  equal  °x  <  °D  '  and  °t  '  resPectivelv 

the  relationship  given  is 


a2   m    (WM.,?  0d2  +  («£da.>»  „t2  (33) 

Since  a    2    can  be  set  equal  to  zero  if  the  measurement  of  t  has 
negligible  error,  the  variance  of  D  calculated  from  N  data  points 
is  given  as 

2      2.    r   ,3f(t,D),2.-l  (34) 
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2 
Brownlee  (9)  has  presented  a  theory  for  an  estimate  of  ox 

based  upon  the  residuals  between  the  least  squares  fit  and  the 

data.   The  best  estimate  of  the  variance  of  X  is  given  by 


x 


2  =  (N-c)"1  I    (Xei-Xei)2  (35) 

i=l 


2 
where  (N-c)  =  number  of  degrees  of  freedom  associated  with  ox 

c  =  number  of  variables,  which  was  one  for  this  analysis. 

Thus,  from  Eqs.  (34)  and  (35)  the  standard  deviation  for  D  is 


r  ?  ,af(t,D),2,-l}l/2   (36) 


where  the  dependence  upon  the  number  of  data  points  N  is  expli- 
citly noted  and  before  the  uncertainty  arising  from  pressure 
measurements  is  considered. 

It  should  be  noted  that  the  measurement  of  pressure  does  not 
enter  directly  into  the  determination  of  D.   However,  pressure  is 
the  independent  variable  reported  with  D  values  and  with  which 
they  are  plotted.   Although  results  could  be  reported  showing 
uncertainties  in  both  the  diffusion  coefficient  and  the  pressure, 
it  is  desirable  to  report  them  in  the  conventional  manner  with  an 
uncertainty  associated  entirely  with  the  dependent  variable,  D. 
To  determine  this  equivalent  description  of  variance,  it  is 
assumed  that  the  total  variance  of  D  is  the  sum  of  that  originating 
in  the  least  squares  analysis  and  the  variance  associated  with  the 
uncertainty  in  the  pressure  measurements.   An  estimate  of  the 
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latter  quantity  is  obtained  by  assuming  Bosanquet's  relationship 
to  be  true 

-1 


D    =    Deil    -     [1/DR      +    1/D-q] 


=  %  db/(dbi  +  DKlp»  (37) 

2 

where  D„   =  the  Knudsen  diffusion  coefficient,  cm  /sec 
Kl 

D_.  =  the  diffusion  coefficient  at  atmospheric  pressure, 
Bl 

2. 
cm  /sec 

P  =  the  pressure,  atm. 

Assuming  a  normal  propagation  of  error,  Eq.  (33)  may  then  be 

applied  to  yield 

°D2(P)  =  [DB1  PKl/(DBl  +  DK1P)2]V 

=  [DeVV  V  •  (38) 

The  value  of  D  , ,  is  obtained  from  Eq.  (37). 
ell 

2 
In  order  to  estimate  the  values  of  op  ,  the  absolute  pressure 

was  expressed  in  terms  of  polynomial  functions  of  the  pressure 

gauge  readings.   The  parameters  of  these  functions  and  the 

variances  of  the  parameters  were  found  by  use  of  POLYFIT  (19) ,  a 

least  squares,  polynomial  fitting  computer  program.   The  order  of 

polynomial  to  be  used  for  the  description  of  a  calibration  curve 

was  selected  on  the  basis  of  the  location  of  a  local  minimum  in 

plots  of  quality  of  fit  number  versus  order  of  fit.   It  should  be 

noted  that  actual  pressure  determinations  were  made  directly  from 
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calibration  curves.   Use  of  the  describing  functions  was  re- 
stricted to  estimation  of  the  variance  cf  the  experimental  diffu- 
sion coefficient.   By  use  of  these  functions  and  again  the  assump- 
tion of  a  normal  propagation  of  error,  an  absolute  pressure  value 
and  an  associated  variance  could  be  found  for  a  given  pressure 
gauge  reading. 

Knowledge  of  c  2  and  De;,,  allow  the  evaluation  of  Eq.  (38) 
to  give  an  estimate  of  the  variance  of  D  caused  by  uncertainty  of 
pressure  measurements.   The  deviation  of  the  effective  diffusion 
coefficient  is  finally  given  as  the  square  root  of  the  sum  of  the 
variance  originating  in  the  least  squares  analysis  of  count  rate- 
time  data  and  the  variance  given  by  Eq.  (3  8)  as 

°D  "  f"D2<N>  +  °D2<P>}1/2  (39) 
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3.0   RESULTS  AND  DISCUSSION 
3.1  General  Discussion 

The  purpose  of  this  study  was  to  verify  and  extend  the  work 
of  Mistier  (5  3)  and  to  refine  and  further  demonstrate  the  feasi- 
bility of  the  experimental  and  analytical  techniques  he  employed. 
A  modified  Ney-Armistead  cell  was  used  with  a  tracer  technique 
which  utilized  a  semiconductor  radiation  detector  within  the  cell. 
Effective  gaseous  self-diffusion  coefficients  were  measured  in 
the  transition  pressure  range.   This  pressure  range  was  chosen 
because  of  the  great  scarcity  of  diffusion  data  for  it.   The  re- 
sults given  in  this  section  illustrate  the  success  of  the  study 
in  achieving  these  goals.   In  particular,  the  C02  data  represents 
the  results  of  improvements  in  experimental  apparatus  and  was 
taken  to  substantiate  Mistier' s  results.   In  addition  data  for 
methane  was  obtained. 

Experimentally  determined  values  of  the  self-diffusion  co- 
efficient for  CO   in  the  transition  pressure  range  at  20°C  are 
presented  in  Table  II  and  Fig.  9.   Values  for  CH4  are  shown  in 
Table  III  and  Fig.  10.   The  results  for  each  gas  have  been  com- 
pared with  theoretical  values  predicted  by  the  relationships  of 
Bosanquet  and  Wheeler,  Eq.  (19)  and  Eq.  (24)  respectively.   Values 
of  3;q  in  the  theoretical  calculations  were  derived  from  the 
theory  of  Chapman  and  Cowling  (12)  at  760  mm  pressure  and  20 °C. 
For  CO- ,  the  theoretical  self-diffusion  coefficient  under  these 
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Table  II.   Experimental  and  theoretical  results  for 
the  self-diffusion  coefficient  of  carbon 
dioxide  in  the  transition  pressure  range. 


p 

Microns 
of  Hq 

D  (cm' 
e 

Experit 

self-di 

718.1 

t 
/sec) 

nental 
.ffusion 

2 
D  (cm  /sec) 

Bosanquet 

2 
D  (cm  /sec) 

Wheeler 

46 

+  45.5 

735 

899 

104 

567.5 

+  24.2 

484 

628 

204 

348.2 

+  14.1 

304 

384 

306 

257.7 

i    9-1 

221 

266 

442 

222.3 

+   6.1 

162 

186 

451 

194.1 

+   6.1 

159 

182 
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Figure  9.   Comparison  between  experimental 
and  theoretical  self-diffusion  coefficients  for  CO~ 
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Table  III.   Experimental  and  theoretical  results 
for  the  self-diffusion  coefficient  of  methane 
in  the  transition  pressure  range. 


p 

D  (cm' 
e 

i 
/sec) 

2 
D  (cm  /sec) 

2 
D  (cm  /sec) 

Microns 

Experimental 

Bosanquet 

Wheeler 

of  Hq 

self-di 

.f fusion 

47 

1184.3 

+ 

35.0 

1290 

1570 

52 

1102.8 

+ 

35.0 

1240 

1520 

53 

1162.6 

+ 

45.4 

1230 

1510 

74 

1018.0 

+ 

21.6 

1050 

1330 

78 

963.6 

+ 

25.7 

1030 

1300 

80 

1028.7 

+ 

24.7 

1010 

1290 

102 

988.0 

+ 

37.0 

885 

1140 

108 

842.1 

+ 

20.6 

856 

1110 

127 

777.4 

+ 

20.3 

774 

1000 

15  0 

70  7.0 

+ 

18.1 

694 

900 

157 

656.7 

+ 

19.6 

673 

871 

175 

659.5 

+ 

17.1 

624 

804 

198 

632.1 

+ 

16.5 

571 

730 

203 

614.3 

+ 

14.4 

561 

715 

251 

547.8 

+ 

14.4 

478 

596 

308 

586.1 

+ 

8.0 

406 

495 

352 

445.6 

+ 

7.4 

364 

436 

396 

382.0 

+ 

8.3 

330 

389 

500 

322.9 

+ 

6.2 

270 

309 
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Figure   10.      Comparison   between  experimental 
and   theoretical    self -diffusion   coefficients    for   CH4 . 
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conditions  was  0.102  cm2/sec. ;  for  CH4  it  was  found  to  be  0.203 
cm2/sec.   The  Knudsen  coefficients  were  found  from  Eq.  (11)  for 
the  .1026d»IiD.  capillary  tuba.  They  were  1284  and  2127  cm2/sec 
for  C02  and  CH4,  respectively.   The  reported  values  of  Dell  were 
calculated  by  Eq.  (25)  using  experimental  values  of  Del2 •    the 
diffusion  coefficient  of  the  tracer  gas  in  the  natural  gas.   For 
C140-,  in  C02  the  correction  factor  was  1.011  while  for  C  H4  in 
CH4  it  was  1.029. 

The  deviation  reported  for  each  experimental  diffusion  co- 
efficient was  the  standard  deviation  which  corresponds  to  67% 
confidence  limits.   This  deviation  is  a  measure  of  the  accuracy 
of  the  least  squares  fit  of  the  data  and  of  the  uncertainty  of 
pressure  measurements.   For  the  C02  data  the  error  represented 
by  these  deviations  ranges  from  6.4%  for  the  measurement  at  46 
microns  to  2.7%  for  the  measurement  at  442  microns.   For  CH4  the 
percentage  deviations  range  generally  from  2-3%  with  only  an 
occasional  instance  as  high  as  5.5%  at  pressures  of  approximately 
50  microns.   These  deviations  do  not  include  the  effect  of  error 
in  measuring  system  parameters. 

The  agreement  of  the  experimental  results  for  C02  with 
theoretical  values  is  shown  in  Fig.  2.   At  the  lowest  pressure 
investigated,  46  microns,  the  experimental  value  is  seen  to  be 
in  excellent  agreement  with  the  theoretical  prediction  by  Bosan- 
quet.   At  intermediate  pressures  they  lie  within  the  band  between 
the  two  theoretical  curves.   With  increasing  pressure  the  experi- 
mental values  exceed  those  predicted  by  Bosanquet  by  increasing 


62 


percentages.   This  general  trend  agrees  with  that  found  by  Mistier. 
It  is  thought  to  be  explained  to  some  degree  by  the  occurrence  of 
surface  diffusion.   This  phenomenon  occurs  by  the  mechanism  of 
molecular  transport  with  an  absorbed  layer  on  the  surface  of  the 
capillary  tube.   Transport  by  surface  diffusion  is  complementary 
to  diffusion  occurring  in  the  gaseous  phase  and  its  effect  would 
be  to  increase  the  observed  effective  diffusion  coefficient. 
Discussions  of  this  phenomenon  are  presented  by  Reynolds  and 
Richley  (64,65),  Clausing  (13),  Hill  (35),  and  others  (11,13,18, 
35,57,68,88) . 

The  agreement  of  experimental  results  for  CH4  with  the 
theoretical  curves  is  shown  in  Fig.  3.   The  same  general  trend  is 
noted  in  the  CH4  data.   At  low  pressures  the  agreement  with  Bosan- 
quet's  relationship  is  rather  good  although  the  points  tend  to 
lie  somewhat  below  the  theoretical  curve.   No  explanation  for 
this  tendency  is  readily  available.   At  intermediate  pressures 
the  measured  values  again  lie  within  the  region  between  the  two 
theoretical  curves;  however,  this  behavior  occurs  at  higher  pres- 
sures than  for  CO,.   Finally  at  the  upper  pressure  region  of  the 
transition  range  the  measured  values  consistently  exceed  the 
values  of  Bosanquet's  relationship.   Again  this  is  thought  to  be 
explained  to  some  extent  by  surface  diffusion. 

The  occurrence  of  surface  diffusion  is  thought  to  account 
for  the  largest  single  source  of  experimental  error.   At  the 
present  time,  knowledge  of  this  phenomenon  is  largely  restricted 
to  very  generalized  theoretical  results  (64)  or  to  studies  of 
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very  specialized  cases.   No  attempt  was  made  to  analyze  its 
quantitative  effects  in  this  study. 

Pressure  control  and  measurement  also  were  responsible  for 
some-  error  in  the  determinations  which  is  represented  by  Oq(P) 
and  is  included  in  the  reported  standard  deviation.   This  uncer- 
tainty in  pressure  measurement  generally  accounted  for  less  than 
30%  of  the  total  deviation.   The  procedure  for  this  analysis  is 
explained  in  Section  2.6.2. 

The  presence  of  a  pressure  gradient  between  the  two  flasks 
would  of  course  result  in  an  erroneous  measurement  of  the  diffu- 
sion coefficient.   By  use  of  pressure  gauges  in  both  flasks  this 
problem  was  essentially  eliminated.   Analysis  of  this  source  of 
error  by  the  approximation  presented  in  Appendix  E  indicated  that 
it  contributed  an  error  of  considerably  less  than  1%  to  the  re- 
ported values. 

Mistier  (53)  has  presented  a  study  of  errors  caused  by  in- 
accuracies in  the  measurement  of  system  parameters  and  in  general 
these  errors  were  found  to  be  small  in  comparison  with  the 
counting  error  associated  with  the  reported  standard  deviation. 

A  final  possible  source  of  error  is  the  Rayleigh  end- 
correction  which  was  used  for  determining  the  effective  tube 
length.   Although  this  is  the  accepted  method  of  correction,  it 
does  not  allow  for  expected  changes  in  the  diffusion  length  as 
Knudsen  diffusion  becomes  predominant.   Presently  no  relationship 
has  been  accepted  for  the  change  in  effective  tube  length  as  a 
function  of  pressure. 


64 


3.2   Conclusion 

Comparison  of  the  results  presented  above  with  the  stated 
purposes  of  the  study  allow  some  conclusions  to  be  drawn. 

Refinements  in  Mistier 's  experimental  equipment  have  re- 
duced the  importance  of  and,  in  some  cases,  eliminated  sources  of 
potential  error.   Refinements  in  the  analysis  of  data  have  per- 
mitted more  meaningful  estimates  of  errors.   The  small  amount  of 
C02  data  presented  was  in  good  agreement  with  Mistier' s.   This 
indicates  that  results  are  reproducible.   The  feasibility  of  using 
a  semiconductor  detector  in  a  modified  Ney-Armistead  diffusion 
cell  to  measure  diffusion  coefficients  has  been  further  demon- 
strated by  the  investigation  of  the  C14H4-CH.   system.   Effective 
self-diffusion  coefficients  of  CH4  at  pressures  from  50  to  300 
microns  were  found  to  compare  within  10%  with  the  theory  of 
Bosanquet. 

Although  the  refinement  of  this  method  must  be  considered 
incomplete,  the  results  of  this  study  suggest  that  further  ef- 
forts will  produce  a  useful  research  tool  for  determination  of 
effective  diffusion  coefficients  in  the  transition  pressure  range. 
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4.0   SUGGESTIONS  FOR  FURTHER  STUDY 

A  first  step  in  the  extension  of  this  study  would  be  a 
thorough  investigation  of  the  effects  of  surface  diffusion.   It 
is  believed  that  modification  of  the  mathematical  model  to  in- 
clude these  effects  would  eliminate  a  major  source  of  error.   The 
use  of  more  precise  pressure  measuring  devices  would  be  highly 
desirable.   The  determination  of  system  parameters  could  be  im- 
proved by  measuring  the  volumes  of  the  flasks  after  final  assembly 
with  a  gas  instead  of  a  liquid;  the  capillary  tube  radius  could 
be  determined  more  accurately  by  optical  or  viscosity  measure- 
ments, particularly  if  its  bore  were  smaller.   Efforts  should  be 
made  to  obtain  more  precise  temperature  control,  preferably  +0.01° 
C.   If  possible,  future  studies  should  employ  tracer  gas  of  a 
higher  specific  activity  to  reduce  errors  produced  by  low  count 
rates.   The  length  of  the  capillary  tube  should  be  varied  in 
order  to  extend  the  pressure  range  to  be  investigated.   This 
would  also  allow  the  determination  of  an  end-correction  for  the 
tube  length  as  a  function  of  pressure. 

Other  possibilities  for  modification  of  the  apparatus  in- 
clude the  installation  of  detectors  in  both  flasks.   This  change 
coupled  with  a  mathematical  solution  in  which  the  tube  volume  is 
considered  negligible  would  greatly  simplify  data  analysis,  and 
allow  several  measurements  within  a  run.   It  would  however  compli- 
cate the  supporting  electronic  system. 

An  exact  mathematical  solution  or  method  of  numerical  anal- 
ysis for  the  system  under  conditions  of  simultaneous  diffusion 
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and  pressure  induced  flow  would  not  only  simplify  the  demands  of 
the  experimental  procedure,  but  would  also  allow  the  measurement 
of  other  transport  properties,  such  as  an  effective  transition 
range  viscosity. 

The  possibilities  for  study  of  other  gaseous  systems  are 
almost  unlimited.   Included  in  these  are  the  measurement  of  self- 
diffusion  coefficients  for  other  gases,  the  measurement  of  mutual 
diffusion  coefficients  for  binary  mixtures,  and  determinations  of 
the  concentration  dependence  of  mutual  diffusion  coefficients. 
By  use  of  multiple  tracers  and  pulse  height  analysis  studies  could 
be  performed  in  multicomponent  systems. 
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7.1   APPENDIX  A 

Solution  to  the  Differential  Equation  for 
Diffusion  Through  a  Capillary  Tube 

Based  upon  the  assumption  that  Fick's  second  law  is  valid 
for  the  process  of  gaseous  diffusion  in  a  capillary  tube,  the 
partial  differential  equation  which  described  gaseous  diffusion 
in  the  system  under  study  was 

3C  _    a2C  ,,  ,, 

o  Z 

In  this  section,  Eq.  (A-l)  is  solved  by  the  method  of  Laplace 
transforms  for  the  three  following  boundary  conditions: 

t=0,   z  >  0,   C=CX  (A-2) 


'0    "  ^fl 

°    1  In 


n0(-D*R2)f§|z=o 


t    >    0,       Z    =    L,       C|z=L    =    C1    +    V-^JVd.R2)     §§|g?J 

nR2Cdz   +   C,irR2L)  (A-4) 

0  L 


To  transform  Eq.  (A-l)  and  the  boundary  conditions  to  dimen- 

2 

sionless  equations  let  X  =  C/C ,  5  =  z/L,  and  t  =  tD/L  .   These 

definitions,  when  applied  to  the  format  of  Eqs.  (A-l) ,  (A-2) , 
(A-3)  ,  and  (A-4)  ,  yield 

ixir^il  =  a!Mu£L    ,  (A_5, 

3T  ar 
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while  the 

boundary  conditions  become 

T  =  0, 

5  >  0,   x  =  xx 

(A- 

6) 

T   > 

0,   5=0, 

T 

..x-.JjilBJJ.* 

>  0,   £-1, 

(A- 

7) 

X( 

t,D 

=  (B  +  1)XX 

-  s((T  3X'I'0)  t  ♦  f'xdc) 
J0   3e       Jo 

(A- 

•8) 

whe  re  a  = 

irl^LV"1  and  B  = 

tS  LVl  . 

Equa 

tion 

(A-5) ,  with 

the  accompanying  boundary 

conditions 

(A-6)  ,  (A 

-7), 

and  (A-8)  ,  : 

Ls  transformed  with  respect 

to  the 

independent  variable  t  ,  thus  letting 

L[X(t 

,5) ]  =  X(s,5) , 

and  using 

Eq. 

(A-6) ,  Eq. 

(A-5)  becomes 

sX(s,5) 

„  _  d2X(s,5) 

1        2 
1     dKd 

(A- 

■9) 

Solving 

this  ordinary  second  order  differential 

equation , 

it 

is  found 

that 

X( 

s,5) 

=  A  (s)  sinh  (X?)  +  B  (s)  cosh  (A?)  + 

s 

(A- 

-10) 

where  X  = 

s1/*. 

To  determine 

the  coefficient  functions  A  (s)  and  B  ( 

s)  it  is 

observed 

that 

the  Laplace 

transforms  of  boundary  conditions 

(A- 

■7) 

and  (A-8) 

are 
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T  >  0, 

K   =  0, 

X(i 

3,0)  =  s"1  +  as"1  dX<='°> 
d? 

(A-ll) 

and 

t  >  0,   5 

=  1, 

X(s 

,1)  =  (B  +  DXjs"1  -  6s"1  ^|- 

,0) 

-  B 

X(s,C)dC  . 

0 

(A-12) 

Taking  the  derivative  o 

f  Eq 

.  (A-10)  and  applying  boundary 

condi- 

tion 

(A-ll)  , 

B  (S) 

=  s 

^(l  -  X.)  +  as"1  dX<  =  '0) 
1            d£ 

(A-13) 

=  S^d 

-  V 

+  as"  [A  (s)  X  cosh  (0)  +  B  (s)  sinh  (0)] 

thus 

B 

(s)  = 

(1  - 

X^s"1  +  aA  (s)  X-1. 

(A-13a) 

Apply 

ing  boundary 

condition 

(A-12)  to  Eq.  (A-10) 

A  ( 

s)  sinh  (XI) 

+  [d 

"  X] 

)s-1  +  aA  (s)  X_1]  cosh  (XI) 

+ 

Xj^s"1  =  (8  + 

DXjs' 

-1 

BA  (s)  X-1  -  S[1{A  (s)  sinh  (X£) 
J0 

+  [(l  -  xr 

Is"1* 

aA  1 

s)  X~  ]  cosh  (X£)  +  X.s"  }d£. 

(A-14) 

After 

performing  the  indicated  integration  and  simplification  of 

Eq.  (A-14) ,  the  coefficient 

function  A  (s)  is  found  to  have  the 

form 

79 


A  (s)  =  (X.  -  1)16  sinh  (A)  +  X  cosh  (X)][(s  +  aB>  X  sinh  (X) 

+  (a  +  B)s  cosh  (X)]"1.  <A-15) 

Substituting  the  expressions  for  A  (s)  and  B  (s)  into  Eq .  (A-10) 
and  simplifying  gives  the  solution  for  X(s,E)  as  a  function  of 
the  variables  s  and  £  as  follows: 

X(s,U  =  X^"1  +  (1  -  Xjs"1  cosh  <X0  +  (X1  -  Us^KS  sinh  (X) 
+  X  cosh  (X)][X  sinh  (XO  +  a  cosh  (XC)1 
•  [(s  +  aB)  sinh  (X)  +  (a  +  B)X  cosh  (X)]"  (A-16) 

In  order  to  use  one  of  Heaviside's  expansion  theorems  in 
finding  the  inverse  of  Eq.  (A-16),  it  is  necessary  to  find  the 
unrepeated  roots  of  both  the  first  and  second  terms  of  the  equa- 
tion.  Clearly,  zero  is  a  root  of  both  terms  and  only  the  second 
term  has  non-zero  roots  of  higher  order.   These  roots  are  found 
by  equating  the  denominator  of  the  second  term  to  zero  as  follows: 

at(o6  +  s)  sinh  X  +  (a  +  B) X  cosh  X]  =  0       (A-17) 

which  leads  to  the  expression 

tanh  Xn  =  -Xn(a  +  8)/(aB  +  X^) .  (A-18) 

Equation  (A-18)  has  an  infinite  number  of  imaginary  roots  desig- 
nated by  Xn. 

Finding  the  inverse  of  Eq.  (A-16)  by  Heaviside's  formula 

using  the  roots  given  above  yields 
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i    i st    .                          t    ~nd    . 

X(t,«)     =    L"1!1         term    1    +    L_1[2         termJ 
zero   root                    zero   root 

+  L-lt3rd   term              -1        3rd   term 

zero   root                    negative   roots 

(A-19) 

where 

T      1  st     . 

-1.1        term   .    _ 
U      'zero   rootJ         xl 

(A-20) 

.    ~nd   . 
-1.2        term   ,         ,         „ 

L      'zero   root'         X        Xl 

(A-21) 

i    -,rd    . 
-1.3        term   ,    _       .    i 

L      'zero   root1    "   0/*|A2   =   Q 

(A-22) 

where 

6   =    (Xx    -    1)  [B    sinh    (A)    +    A    cosh    (A)  ] 

2 
•     [A    sinh    (AS)    +    a    cosh    (A£)]e      T 

*    =    (oB    +    A2)     sinh    (A)    +    A(a   +    8)    cosh    (A)    +    A2 

-1                    2 
•     [(2A)        (oB    +    A      +    a    +    B)    cosh    (A) 

+    (l/2(a   +    6)    +    1)    sinh    (A) ] 

(A-23) 

so   that 

0/<fi|    ,         =    a(X,    -    1)  (B    +    1)  (a   +    B    +    aB)"1 

|a  =0             l 

L      'negative'roots'    =    <X1   "    «    I,    en/*n 
^                                                      n=l 

(A-24) 

where 
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0   =  [B  sinh  Un)  +  Xn   cosh  Un>] 
■     [X   sinh  (X  £)  +  a  cosh  (Xn£)]e     , 


i*   =  X2[(2X  )-1(a  +  B  +  aB  +  X2)  cosh  (X  ) 
n    n    n  n         n 

+  (1/2  (a  +  8)  +  1)  sinh  (Xn>  ] 

2 
Equation  (A-24)  is  simplified  by  redefining  Xfl  to  have  a 

positive  value.   This  operation  requires  the  substitution  of 

-X2  for  X2  in  Eqs.  (A-18)  and  (A-24).   Using  the  identities 

n      n     ' 

tanh  ix  =  i  tan  x, 
cosh  ix  =  cos  x, 
sinh  ix  =  i  sin  x, 

and  the  definition  for  Xn ,  causing  Xn  to  equal  iXn,  Eqs.  (A-18) 
and  (A-24)  reduce  to 

tan  Xn  =  -Xn(a  +  B)/(a8  -  X2)  (A-25) 

and 

-1    3 r ^  term    ,  =  _  f    t (      , /  t (  (A_26) 

'negative  roots      *,   n  n'      Tn  n' 

where 

0T(X  ,£,t)  =  e   n  [a  cos(XE)  -  Xn  sin  (X  £)] 
n   n  n       n       n 

•  [X   cos  X   +  6  sin  X  1 
n      n  n 


82 


*n(Xn>      =   Xn{(2Xn>       (a   +    B   +   aB    "    Xn>    cos    Xn 
-     [1/2  (a   +    8)    +    1]    sin    Bn ) . 


n 


The  eiqenvalues  X   of  Eq.  (A-25)  are  solved  accurately  by  Newton's 

3  n 

method  (33)  given  a  close  approximation  to  the  actual  Xn  values. 
The  first  order  solution  X-^   is  given  by 

tan  X,  «  X,  -  X1(a  +  6)/.Uj  -  aB)  (A-27) 

1/2 

so  that  X,  «  (a  +  8  +  aB)  '  . 

Approximations  to  higher  order  values  of  Xn  are  obtained  by  use 


of  the  algorithm  X   *  xn-l  +  *" 

Thus  the  solution  to  the  partial  differential  equation  (A-5) , 
given  by  Eq .  (A-19) ,  can  be  written 

X(t,5)  =  (B  +  aX^B  +  l)/(a  +  B  +  08) 

+  fXr  "  1)  I     ♦B(Xn'«'T,/*n{xn)  (A_28) 

n=l 

where  X   is  a  solution  to  the  transcendental  equation  (A-25) . 
n 

Equation  (A-27)  shows  the  change  in  concentration  gradient  at  a 
distance  variable  point  z  with  respect  to  time  as  a  function  of 
the  system  parameters  V,,  V2 ,  R,  L,  and  D. 
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7.2   APPENDIX  B 

Leakage  Error  Analysis 

In  order  to  correct  for  leakage  of  the  tracer  gas  prior  to 
opening  of  the  capillary  valve,  "10",  Mistier  subtracted  a  "back- 
ground" value  from  all  experimental  count  rates  in  each  run.   This 
subtraction  forced  the  normalized  count  rates  to  fit  the  boundary 
conditions  of  his  mathematical  model.   The  model  assumed  zero 
concentration  of  the  tracer  gas  in  the  capillary  tube  and  right- 
hand  flask  at  the  time  of  opening  of  the  capillary  valve.   The 
"background"  count  rate,  k.,  was  obtained  by  solving  the  equation 

(C   -  k.)(C   -  k,)"1  =  6(a  +  e  +  aB)"1  (B-l) 

co      x     o      X 

where  C  =  the  observed  equilibrium  count  rate, 

C  =  the  observed  initial  count  rate, 
o 

a  and  8  are  defined  as  in  Appendix  A, 

8(a+B+a8)    =  the  normalized  equilibrium  concentration  for 
an  initial  zero  concentration  of  tracer  gas  in  the  tube  and 
right-hand  flask. 
If  a  quantity  Xm  is  defined  as 

X°  i  Bta+B  +  aB)"1  <B"2) 

the  solution  of  Eq.  (B-l)  for  k.  yields 

k,  =  (C   -  C   X°)(l  -  Xf)"1  (B-3) 
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If  it 

is  noted  that  in  most  circumstances  only  the  first  term  of 

the  summation  makes  a  significant  contribution,  an  approximate 

diffui 

sion  coefficient  is  found  to  be 

D'  =  -L2(A^ti)"1  lM-fX^-X^A"1 

(B-7) 

If  in 

the  solution  as  found  in  Appendix  A  the  same  approximation 

of  ut: 

Llizing  only  the  first  term  of  the  series  is  applied, 

the 

diffusion  coefficient  becomes 

D  =  -L2(s,t.)_:L  ln(-(X.-X  )  (1-XC)  (A.  (1-X  ))_1) 

1  l           i   ■      "    1     " 

(B-8) 

Based 

on  these  assumptions,  the  error  associated  with  Mistier 's 

approximation  is  for  the  i   observation 

D-D 

'  =  -L2(s,t.)  1   InU.-X  )(l-xC)(d-X  )  (X.  -XC))_1) 

li        i°°     ™      "   la   " 

(B-9) 

The  computer  code  ERCALC  was  used  to  compute  the  value  of  this 

error 

As  shown  in  Appendix  A,  the  normalized  count  rate. 

KL,    at 

time, 

t . ,  is  given  by 

X,  =  X  +  (X,-l)  I      A  e 
1           l           n=l   n 

(B-10) 

where 

X.  and  t.  are  defined  as  in  Appendix  A.   It  is  noted 

that 

t.  is 

l 

contained  in  i.  as 

xx;  t.  dl-2. 

(B-ll) 

85 

Performance  of  the  described  subtraction 

gives 

the  approximate 

value 

of 

the  i   normalized  concentration  as 

Xia  = 

(C.  (1-XC)  -  C   +  C  XC)(C  (1-XC)  -  C   +  C 

loo           00         Q   00       Q         00           ooc 

X0)"1 

»  oo' 

IB- 

-4) 

in  terms 

of  concentrations  normalized  by 

the  observed 

initial 

count 

rate  this  is 

Xja  =  'V1-*"'  _  Xoo  +  0<1_X. 

y1 

1 

(B- 

■5) 

where 

X. 

1 

is  the  i   normalized  concentration  as 

used  : 

Ln  the  so] 

.u- 

tion  : 

tor 

a  finite  initial  concentration  in  the 

capillary  tube  and 

right- 

-hand  flask  (Appendix  A)  . 

Mistier' s  model  was 

«      -t. 

X,  =  =  B(a+B+aB)    -   7   A„e 
la                 L,   n 
n=l 

A^D/L2 

• 

(B- 

■6) 

where 

X2t 
An  =  [8nUn'C'T)/*n(V^  "T 

and 

»^V£|>  »♦„<»»>  ' 

and 

X      are  defined  in  Appendix  A. 
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Representative  values  of  X^  and  X.  were  chosen  and  used  as  input 

2 

values  to  ERCALC.   Values  of  D  ranging  from  100  to  1000  cm  /sec 

were  employed  as  were  values  of  t.  for  each  value  of  D.   The 
values  of  t.  were  chosen  so  as  to  span  the  transient  time  of  the 
system  for  that  D. 

The  subroutine  XCALF  was  called  within  ERCALC  for  each  D  and 
its  associated  t.'s,  and  by  use  of  Eq.  (B-10)  values  of  X^^  were 
generated.   ERCALC  then  computed  the  error  of  Mistier' s  approxima- 
tion by  Eq.  (B-9)  at  each  value  of  ti  for  each  D. 

The  maximum  magnitude  of  this  error  in  all  cases  investigated 
was  found  to  be  considerably  less  than  1%.   This  indicated  that 
the  error  induced  by  Mistier' s  approximation  was  negligible  when 
compared  with  the  other  sources  of  error  he  encountered. 
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7.3   APPENDIX  C 

Data  Analysis  for  the  Diffusion  Coefficient 

7.3.1  Computer  Program 

A  computer  program  previously  developed  by  Mistier  for  anal- 
ysis of  data  from  this  system  was  modified  to  agree  with  the 
solution  for  the  system  with  an  initial  concentration  of  tracer 
gas  in  the  capillary  tube  and  right  hand  flask  (See  Appendix  A) . 
Basically,  as  in  its  original  form,  the  program  performs  an 
iterative  least  squares  analysis  to  determine  the  "best"  value 
for  the  effective  diffusion  coefficient.   Minor  additions  were 
made  to  allow  internal  normalization  of  data  points  and  the 
association  of  a  time  with  each  point.   A  subroutine  was  added  to 
the  program  to  calculate  the  increase  of  the  variance  of  the  ex- 
perimental diffusion  coefficient  caused  by  uncertainty  in  pressure 
measurements  (See  Section  2.6.2).   A  brief  discussion  of  the 
least  squares  analysis  follows. 

7.3.2   Least  Squares  Analysis  and  Numerical  Approximations 

In  the  method  of  least  squares  (9)  the  sum  of  the  squares  of 
deviations,  Re,  between  the  observed  values,  Xei(  and  the  Xc. ,  are 
minimized  to  find  the  values  of  Xc.  which  give  the  minimum  error. 
Applying  this  method  of  analysis  to  the  diffusion  data  the  ex- 
pression for  the  residue,  Re,  is 
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Re  =  I    (Xe,  -  Xc.)' 
i=l   1     1 


(C-l) 


where  Xe .  and  Xc.  are  the  ith  experimental  and  calculated  values 
of  X  at  time  t.  respectively,  and  N  is  the  number  of  data  points. 
The  minimum  squared  error  is  found  by  differentiating  Re  with  re- 
spect to  D  and  setting  the  resultant  equal  to  zero;  thus, 


3  Re 
3D 


N  3Xc. 

2  I     <*4  -  XC,)  -jgi-  0 

1  =  1 


(C-2) 


Defining  a  quantity 


N  3Xc. 

fD=  .1  J^-Xc.)  -jgi 


i=l 


and  utilizing  the  first  two  terms  of  Taylor's  expansion  (33) 
the  result  is 


fD  "  fD   +  "  15° 
o 


(C-3) 


where  f   should  be  zero  by  Eq.  (C-2)  so  that 


—  3fD   _1 

o 


(C-4) 


To  evaluate  the  right  hand  member  of  Eq.  (C-4)  use 


3D 


=  I 
i=l 


(Xe,  -  Xc,) 


32Xc.       3Xc. 

1  _  /  ±_   \ 

3D^         3D   ' 


(C-5) 


where 


3Xc-     N   2  /T2.  t ,    .   ,  „t 


30    i=l 


32Xc, 
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I  (<tA')»;(»n,5.t)/*;(v  <c-6> 


r"    J    (^TA2)29>n,i,T)A>n)  .  (C-7) 


3D        i=l 

where  X  ,  x,  L,  I,    6T,  and  $f   are  defined  as  in  Appendix  A.   The 
value  of  D  with  minimum  error  can  be  determined  by  using  the 
algorithm  D   ,  =  D  +  AD  and  iterating  upon  D  until  AD  is  suffi- 
ciently small. 
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7.4   APPENDIX  D 
Experimental  Procedure 

Prior  to  the  actual  accumulation  of  data,  several  prepara- 
tions were  necessary.  After  carefully  installing  the  surface 
barrier  detector,  all  valves  except  "15",  Fig.  (5),  were  left 
open  with  the  vacuum  pump  running  for  approximately  2  4  hours. 
This  allowed  complete  desorption  of  water  vapor  adsorbed  on  the 
glass  surfaces  of  the  system  during  exposure  to  the  atmosphere 
and  attainment  of  a  pressure  of  approximately  ly  of  Hg.   After 
degassing,  the  rate  of  desorption  was  so  slow  that  there  was  no 
detectable  pressure  change  during  the  time  required  for  an  ex- 
periment.  This  time  ranged  from  20  to  60  minutes  approximately. 

The  system  was  then  purged  with  CH4  or  C02  a  minimum  of  six 
times.   For  each  purge,  the  pressure  throughout  the  system  was 
raised  to  approximately  500  mm  of  Hg  by  closing  valve  "9"  and 
bleeding  the  appropriate  gas  in  through  valve  "15".  Valve  "15" 
was  then  closed,  the  vacuum  pump  reintroduced,  and  the  system 
evacuated  to  ly  of  Hg. 

After  purging,  valves  "2",  "3",  and  "4"  were  closed  and  the 
radioactive  gas  flask  removed.   The  glass  stem  of  the  ampoule  was 
broken  with  a  short  glass  rod  housed  in  the  charging  flask  and 
the  flask  replaced.   Valves  "3",  and  "4"  were  opened  and  the  sys- 
tem again  thoroughly  purged. 
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To  calibrate  the  thermocouple  vacuum  gauges,  the  pressure 
was  again  raised  to  approximately  500  mm  of  Hg  by  introducing  CH4 
or  CO, ,  and  the  system  evacuated  to  the  upper  portion  of  the  pres- 
sure range  of  interest,  approximately  500 p  of  Hg.   After  a  period 
of  approximately  three  minutes  (valve  "11",  "5"  and  "12"  remained 
open  throughout  the  calibration  to  afford  rapid  pressure  equaliza- 
tion) to  allow  establishment  of  pressure  equilibrium  throughout 
the  system,  pressure  readings  were  recorded  from  the  electrical 
vacuum  gauges,  "6"  and  "18",  and  from  the  McLeod  vacuum  gauge. 
The  pressure  of  the  system  was  lowered  a  few  microns  of  Hg  and 
after  approximately  three  minutes  the  gauge  readings  were  recorded. 
These  steps  were  repeated  until  the  pressure  was  approximately  one 
micron  of  Hg.   The  pressure  was  again  raised  to  approximately  500 
mm  of  Hg  and  another  traverse  of  the  pressure  range  made.   This 
procedure  was  continued  until  enough  readings  had  been  amassed  to 
construct  reliable  calibration  curves  for  the  thermocouple  vacuum 
gauges . 

The  integrity  of  the  vacuum  was  checked  daily  and  occasional 
spot  checks  were  made  of  the  calibration  curves.   The  thermocou- 
ple gauges  showed  no  drift  from  day  to  day,  but  the  reproducibi- 
lity of  the  readings  varied  somewhat  during  alternate  charging 
and  evacuation  of  the  system. 

In  order  to  guard  against  drift  in  the  electronic  equipment, 
all  components  except  the  oscilloscope  were  supplied  power  by  a 
Beckman  line  voltage  regulator.   Ground  loops  were  prevented  by 
connecting  all  plug  ground  lines  to  receptacle  grounds  and  by 
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connecting  all  chassis  to  a  common  ground.   To  prevent  the  detec- 
tion of  ambient  electromagnetic  noise  signals  flask  I  was  enclosed 
in  a  wire  mesh  which  was  electrically  grounded. 

To  determine  usable  settings  for  the  electronic  equipment, 
valves  "11"  and  "9"  were  closed  and  flask  I  charged  to  a  pressure 
of  approximately  75 p  of  Hg  with  radioactive  gas.   With  arbitrary 
settings  of  the  amplifier  gain  and  postamplif ier  gain  and  a  low 
postamplifier  bias  setting  the  output  signal  was  monitored  with 
the  oscilloscope.   The  detector  bias  voltage  was  then  raised  in 
small  increments  and  the  amount  of  noise  in  the  output  signal  ob- 
served on  the  oscilloscope.   After  an  initial  increase,  the  noise 
level  decreased  to  a  minimum  at  an  indicated  detector  bias  of  10  0 
volts.   Combinations  of  settings  for  the  amplifier,  postamplifier, 
and  multiscaler  gains  and  the  postamplifier  bias  were  investigated. 
Settings  of  x8  for  the  amplifier  gain,  x4  for  the  postamplifier 
gain,  x3.0  for  the  multiscaler,  and  1.10  for  the  postamplifier 
were  selected  on  the  basis  of  maximum  count  rate  with  minimum 
noise.   The  frequency  of  the  square  wave  generator  signal  was 
calibrated  with  a  stop  watch. 

The  actual  diffusion  data  was  taken  by  the  following  steps: 

1.  A  time  increment  for  the  square  wave  generator  suitable 
for  the  desired  pressure  was  chosen  and  set  on  the 
generator  by  stopwatch. 

2.  The  water  level  in  the  constant  temperature  bath  was 
lowered  below  the  level  of  the  capillary  valve.  The 
flap  in  the  wire  screen  shield  was  opened  and  the  tape 
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on  the  flask  removed.   The  capillary  valve,  "10",  «as  then 
closed  and  visibly  inspected.   The  tape  was  replaced, 
the  screen  flap  closed,  and  the  bath  refilled. 

3.  The  system  was  evacuated  to  a  pressure  of  approximately 
lp  of  Hg  and  valve  "9"  closed. 

4.  Valve  "15"  was  used  to  bleed  CH4  or  C02  into  the  system 
and  raise  the  pressure  to  a  level  approximately  75p  of 
Hg  less  than  that  desired  for  the  experiment. 

5.  Valves  "11"  and  "12"  were  closed,  and  gas  bled  into  II 
until  its  pressure  reached  that  desired  for  the  deter- 
mination. 

6.  Valves  "2"  and  "4"  were  opened  and  radioactive  gas  bled 
into  flask  I  through  needle  valve  "3"  to  bring  it  to  the 
desired  pressure.   Valves  "2",  "3",  and  "4"  were  then 
closed. 

7.  The  detector  voltage  was  raised  to  lOOv.  and  the  oscillo- 
scope turned  on. 

8 .  After  a  period  of  two  to  three  minutes  the  readings  of 
gauges  "6"  and  "18"  were  recorded  and  "5"  and  "13"  closed. 

9.  The  output  signal  was  monitored  on  the  oscilloscope  as  a 
check  on  the  signal  quality. 

10.  The  multiscaler  was  started  and  at  the  time  it  entered 
channel  number  11,  the  capillary  tube  valve  was  opened. 

11.  After  the  multiscaler  terminated  counting,  the  trace  was 
printed  out  with  the  automatic  printer. 
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12.  The  trace  was  visibly  inspected  on  the  multiscaler  for 
reliability  and  attainment  of  equilibrium. 

13.  If  equilibrium  was  not  reached  in  the  initial  trace,  a 
second  trace  was  made  to  record  the  equilibrium  value  and 
printed  out. 

14.  Valves  "5"  and  "13"  were  opened  and  after  two  to  three 
minutes  the  final  pressure  readings  were  recorded. 

15.  The  system  was  then  purged  in  preparation  for  the  next 
determination . 
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7.5  APPENDIX  E 

Error  Associated  with  Pressure  Gradients 

Although  efforts  were  made  at  the  start  of  each  experimental 
determination  to  insure  constant  pressure  throughout  the  system, 
there  is  reasonable  cause  to  believe  that  this  was  not  always 
accomplished.   Small  differences  in  the  initial  and  final  readings 
of  the  pressure  gauges  occassionally  occurred.   Even  in  the  ab- 
sence of  these  changes  in  pressure  indications,  the  small  scatter 
in  the  calibration  curves  would  lead  to  reasonable  doubt  of  the 
assumption  of  constant  initial  pressure  throughout  the  system  in 
each  determination.   The  presence  of  an  initial  pressure  difference 
between  the  two  flasks  would  cause  a  transfer  of  material  in  ad- 
dition to  that  resulting  from  diffusive  flow  along  the  concentra- 
tion gradient.   The  direction  of  this  pressure-induced  flow  may 
be  in  the  same  or  opposite  direction  as  the  diffusive  flow  depending 
upon  the  relationship  of  the  total  pressure  gradient  and  the  con- 
centration gradient  of  the  diffusing  material.   A  study  of  super- 
imposed pressure-induced  and  diffusive  flow  in  graphite  has  been 
presented  by  Evans,  Truitt,  and  Watson  (23);  an  analysis  for  a 
capillary  tube  is  given  by  Otani,  Wakao,  and  Smith  (56) .   Based 
on  earlier  work  by  Evans,  Watson,  and  Mason  (22),  Rothfeld  (66), 
and  Scott  and  Dullien  (67)  ,  they  state  the  total  flux  of  component 
1  in  a  binary  system  of  1  and  2  in  the  presence  of  both  a  total 
pressure  gradient  and  a  gradient  in  the  partial  pressure  of  1  as 


96 


N* 


dCl  ,-1  ,  -   ,«  D12, 


S  -  "D12  d^  t1-'1-™'^  +  Dl2/DKll    +  FVl(1+  D^1   (E-1} 


where  D,?  =  the  mutual  diffusion  coefficient  according  to  bulk 

2 

diffusion  theory,  cm  /sec. 

C,  =  the  concentration  of  1,  g  -  mole/cm 
z  =  distance  in  direction  of  diffusion,  cm 
m  =  the  square  root  of  the  ratio  of  the  molecular  weight 
of  1  to  that  of  2 

y.  =  the  mole  fraction  of  1 

2 

D„,  =  the  Knudsen  diffusion  coefficient  of  1,  cm  /sec 
Kl 

2 
F  =  flux  caused  by  a  pressure  gradient,  g-mole/sec-cm 

The  value  of  F  is  given  to  be 

9  dP 

F  =  -[r  P/8p+4rRT/(3(M1v1y1+M2v2y2)) ](1/RT)  §j  (E-2) 

where  r  =  radius  of  the  capillary  tube,  cm 

2 
p  =  total  average  pressure,  dynes/cm 

y  =  gas  viscosity,  poises 

R  =  gas  constant,  ergs/g-mole-°K 

T  =  temperature ,  °K 

MjM.  =  molecular  weights  of  1  and  2,  respectively,  gm/g-mole 

v.;v2  =  mean  molecular  speeds  of  1  and  2,  cm/sec 

y2  =  mole  fraction  of  2 . 

In  the  case  of  isotopic  self-diffusion  m  is  approximately  unity 

and  Eq.  (E-l)  reduces  approximately  to 

dC, 

NI  =  "Dell  -df+    (Dell/Dll)Fyi  (E"3) 
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where  D-,,  =  effective  diffusion  coefficient  according  to  Bosanquet, 

2  , 
cm  /sec 

D,   =  self  diffusion  coefficient  according  to  bulk  diffusion 

2 
theory,  cm  /sec 

1  =  subscript  denoting  tracer  gas. 

Under  these  conditions  F  becomes 

F  =  -(r2P/8p+4rRT/3Mv)  (1/RT)  (AP/L)  (E-4) 

where  L  =  capillary  tube  length,  cm. 

The  assumption  of  a  linear  pressure  gradient  in  the  tube  allows 


replacement  of  ||  by  AP/L,  where  AP  is  the  pressure  difference  of 


the  two  flasks  in  the  system  used  in  this  study  UP  >  0  when  the 
pressure  in  flask  I  is  the  greater) . 

By  the  definition  of  F  and  AP  it  may  be  stated  that 

^|£i  -  (Deil/D11)ART(l/V1+l/V2)F  (E-5) 

2 
where  A  =  cross  sectional  area  of  capillary  tube,  cm 

and  V,,  V3  refer  to  the  volumes  of  the  flasks. 

By  combining  Eqs .  (E-4)  and  (E-5) 

diA£)  =  -k,AP  (E-6) 

dt        1 

ere  k,  =  (Deil/Dl:L)  (rP/8p+4rRT/3Mv)  (A/L)  (l/V^l/Vj)  .   The  solu- 


wh 

tion  of  Eq.  (E-6)  with  appropriate  boundary  conditions  gives 

-k,t 


iP(t)  =  AP0  e 


1  (E-7) 
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In  terms  of  an  effective  half-life  then 

Tl/2  =  °-693/kl    '  (E_7A) 

The  quantity  T,-  is  by  definition  the  time  required  for  an  ini- 
tial pressure  difference  to  decrease  to  one  half  its  original 
value.   Experimental  evaluation  of  T1/2  by  the  above  equations 
shows  that  it  ranges  in  value  from  approximately  330  sec.  at  a 
aressure  of  50u  to  approximately  70  sec.  at  500y. 

The  observed  change  in  the  concentration  of  tagged  material 
in  flask  I  is  directly  proportional  to  the  change  in  the  normalized 
count  rate  which  may  be  represented  as 

aX?  =   AXd  +  axp  <E"8) 

where  AX  =  the  total  observed  change  in  normalized  count  rate, 
ax  =  the  contribution  of  diffusion  to  this  change 


ad 


ax 


'F 


the  contribution  of  pressure-induced  flow. 


The  evaluation  of  these  quantities  requires  the  manipulation  of 
Eq.  (E-3)  to  extract  expressions  for  diffusive  and  pressure- 
induced  flow  as  functions  of  time.   To  avoid  this  procedure  which 
is  quite  complicated  and  of  an  iterative  nature,  it  was  assumed 
that 

AX   =  k  N*   e-1^  (E"9> 

axp  -  K2  oF 

where  k.,  is  a  proportionality  constant  and  Eq.  (E-7)  was  implied. 

Here  N*   is  the  initial  flux  of  tagged  material  leaving  flask  I 

oF 

because  of  the  pressure  gradient. 
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If  it  is  assumed  that  for  a  short  period  of  time  the  contri- 
butions of  diffusive  and  pressure-induced  flow  are  of  the  same 
relative  importance  as  at  time  zero 

AXd  1     m   _od  (E-10) 

AXF  N* 

aA"  lt=0     OF 

where  N*,  is  the  initial  flux  of  tagged  material  leaving  flask  I 
by  diffusive  flow.   Knowledge  of  the  initial  concentrations  of 
tracer  in  the  flasks,  the  assumption  of  a  constant  concentration 
gradient  along  the  length  of  the  tube,  and  an  estimate  of  Den 
allows  the  evaluation  of  N*   by 

dC.  | 
N*   =  _n     ±-  (E-ll) 

od       um  dx    |t=0 

If  the  pressure  and  AP   are  also  known  N*   is  obtained  from 

NoF=  (Dell/Dll'  F*l  IfO    '  (E"12) 

By  use  of  the  system  solution  with  the  estimate  of  D?11  as 
given  in  Appendix  A  the  value  of  AX,  may  be  found  when  t  =  At  and 
Eq.  (E-10)  may  be  solved  for  the  value  of  AXp.   Substitution  of 
AX  and  At  into  Eq.  (E-9)  allows  solution  for  k2 .   All  the  quan- 
tities necessary  to  evaluate  AX  by  Eq.  (E-9)  are  then  known,  and 
the  effect  of  pressure-induced  flow  may  be  removed  from  data  by 
continued  use  of  Eq.  (E-9)  and 


AXd  =  AX^,  -  AXp  (E-13) 
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where  AX_  is  the  observed  change  in  the  normalized  count  rate  at 
time  t . 

The  significance  of  this  correction  was  investigated  by  apply- 
ing it  to  representative  sets  of  data  which  had  already  been  anal- 
yzed.  A  computer  code  was  written  to  do  this  and  sets  of  data 
covering  the  pressure  range  analyzed.   An  estimate  of  APQ  was 
determined  from  the  initial  readings  of  the  pressure  gauges  and 
the  gauge  calibration  curves  and  an  estimate  of  Deil  from  the 
original  least  squares  analysis.   The  original  data  was  corrected 
by  use  of  Eq.  (E-13)  and  again  subjected  to  the  least  squares 
analysis.   In  all  the  cases  investigated,  the  corrected  diffusion 
coefficient  differed  from  that  obtained  from  the  uncorrected  data 

by  less  than  1%.   For  the  sample  data  set  presented  in  Table  I,  this 

2 
correction  increased  the  resultant  Deil  from  614.3  to  617.4  cm  /sec 

an  increase  of  approximately  0.5%. 


ISOTOPIC  GASEOUS  SELF-DIFFUSION 
COEFFICIENTS  AT  TRANSITIONAL  PRESSURES 


by 

GLEN  ROY  CORRELL 
B.  S.,  University  of  Missouri  at  Rolla,  1965 


AN  ABSTRACT  OF  A  MASTER'S  THESIS 

submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree 

MASTER  OF  SCIENCE 

Department  of  Nuclear  Engineering 


KANSAS  STATE  UNIVERSITY 
Manhattan ,  Kansas 


1967 


ABSTRACT 

A  modified  Ney-Armistead  diffusion  cell  utilizing  a  semi- 
conductor radiation  detector  was  used  to  measure  self-diffusion 
coefficients  of  CO,  and  CH4  at  20°C  in  the  transition  pressure 
range.   Summaries  of  existing  theories  of  gaseous  diffusion  in 
the  bulk,  Knudsen,  and  transition  region  were  presented.   The 
existing  methods  of  diffusion  coefficient  measurement  were  summa- 
rized.  The  properties,  characteristics,  and  advantages  of  semi- 
conductor radiation  detectors  were  discussed.   Descriptions  of 
the  experimental  arrangements  and  techniques  and  the  methods  of 
data  analysis  were  given  detailed  treatment.   The  significance 
of  results  was  discussed  and  an  extensive  error  analysis  presented. 

The  diffusion  cell  consisted  of  two  flasks  of  579  ml  and 
512  ml  joined  by  a  capillary  tube  10.00  cm  long  and  0.126  cm  I.D. 
A  surface-barrier  radiation  detector  was  contained  within  the 
large  flask.   The  effective  tube  length  was  determined  by  use  of 
the  Fayleigh  end-correction.1  A  solution  of  the  one-dimensional 
diffusion  partial  differential  equation  was  obtained  by  the  appli- 
cation of  Fick's  second  law  and  the  use  of  Laplace  transforms. 
Data  consisting  of  concentration  versus  time  was  subjected  to  an 
iterative  least  squares  analysis  to  determine  the  effective 
diffusion  coefficient. 

The  semiconductor  radiation  detector  was  found  to  have  the 
advantages  of  internal  detection  of  radioactive  tracers,  rapid 
time  response,  and  high  counting  efficiency.   Refinements  of  the 


apparatus  and  analysis  were  employed  to  reduce  the  errors 

2 
encountered  by  Mistier   and  to  increase  the  value  of  this 

approach  for  measurement  of  effective  diffusion  coefficients  in 
the  transition  pressure  range. 

Values  of  the  self-diffusion  coefficient  of  C02  were  pre- 
sented which  supported  Mistier' s  findings.   Effective  diffusion 
coefficients  were  reported  for  CH.  over  a  pressure  range  of  50 
to  500  microns.   Those  between  50  and  300  microns  agreed  within 
10%  with  the  theory  of  Bosanquet  while  those  from  300  to  500 
microns  exceeded  predicted  values  by  as  much  as  12.3%.   Errors 
were  attributed  primarily  to  surface  diffusion  effects  not  con- 
sidered in  the  analysis.   Standard  deviations  based  on  the 
residues  of  the  least  squares  analysis  and  uncertainties  of 
pressure  measurement  were  generally  2-3% ,  with  isolated  values 
as  high  as  5.5%  near  50  microns. 


1Rayleigh,  Lord,  Theory  of  Sound.   Dover  Publications,  2,  291 
(New  York,  1945. 

2Mistler,  T.  E.,  Radioactive  Tracer  Determinations  of  Gaseous 
Diffusion  Coefficients.   M.S.  Thesis,  Kansas  State  University, 
Manhattan,  Kansas,  1966. 


